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INTRODUCTION 


This report serves as the Pinal Report for research supported by the 
National Aeronautics and Space Administration, Goddard Laboratory for Atmospheric 
Sciences under Grant NSG-5252. (Optimum Employment of Satellite Indirect Soundings 
as Numerical Kodel Input.) The fundamental goal of this project was to Identify 
characteristics of satellite-derived temperature soundings that would signifi- 
cantly affect their use as input for numerical weather prediction models. 

Many of the previous sounding evaluations Involved model impact studies 
that mixed satellite soundings with conventional data, before the error charac- 
teristics of the satellite soundings were fully defined. In contrast, our work 
has emphasized independent evaluations of satellite soundings to better define 
these error characteristics. The article by Koehler presents a Nimbus-6 sounding 
study from February 1976 (during the winter Data Systems Test). Hie rcsulto 
reveal an underestimation of the strength of synoptic scale troughs and ridges, 
and associated gradients in isobaric height and temperature fields. The most 
significant errors occur near the earth's surface and the tropopouee. The instru- 
ments carried aboard Nimbus-6 were prototypes for those carried aboard the TIROS-N 
and NOAA-6 satellites, which provided the global distribution of satellite-derived 
temperature profiles during the important FGGE year observational study (from 
December 1978 through November 1979). 

Soundings from the TIROS-N and NOAA-6 satellites are evaluated in Schmidt, 
Koehler and Horn. Their results are remarkably similar to those from Koehler, 
again showing an underestimation of upper level trough amplitudes leading to 
weaker thermal gradient depictions in satellite-only fields. These errors ahow a 
definite correlation to the synoptic flow patterns. 

Derber, Koehler and Horn used ore of Schmidt et al.'s satellite-only analy- 
sis to initialize a numerical model forecast, and found that these synoptically 
correlated errors are retained in the forecast sequence. If the sounding errors 
were not retained or were more random in nature, they could be more easily com- 
bined with conventional data in data rich regions. However, the correlated nature 
of these sounding errors complicates their incorporation into conventional data 
sets. Thus, this knowledge of the nature of the satellite errors should be used 
in additional efforts to develop netho • . which can faclllitate the Inclusion of 
satellite-sounding Information into Ve observational data mix. 
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A Case Study of Height and Temperature Analyses 
Derived from Nlmbus-6 Satellite Soundings 

Thomas L. Koehler 

Department of Meteorology • 

University of Wisconsin-Madison, WI 53706 


Abstract 

Height and temperature analyses were constructed on a subset of the 
LFM grid using only Nimbus-6 satellite temperature profiles from approximate- 
ly 1800 GMT 22 Feb, 1976. Several experiments were performed to evaluate 
various features of these satellite-derived analyses. Fields derived from 
the bracketing LFM analyses provide the verification data. Results from this 
study provide insight into possible reasons for the inconclusive results from 
the DST-6 impact studies. 

The results indicate that Nimbus-6 soundings were able to correctly 
position the major troughs and ridges, but underestimate gradients in the 
analyses due primarily to the soundings being too warm in the troughs. No 
advantage could be found in using a set of satellite soundings with a greater 
horizontal resolution than the DST soundings. 

While the satellite soundings for this period were 'degraded due to a 
loss of one set of infrared channels, results from this study are quite sim- 
ilar to those presented in Phillips et al. (1979), Schlatter (1981) and Schmidt 
et al. (1981) with more recent TIROS-N and NOAA-6 sounding data. This sug- 
gests inherent limitations in the methods used to derive these soundings, and 
in their incorporation into conventional data sets. 
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1. Introduction 

The ability to provide accurate synoptic scale weather predictions 
from one to several days has been hampered by the lack of conventional data 
over large regions of the globe, especially ocean regions. During the last 
decade considerable effort has been expended to alleviate this problem by 
providing truly global data sets. The culmination of this effort was reached 
in the collection of the FGGE year data sets (December 1978 to November 1979). 
Vertical temperature profiles derived from satellite measurements comprise an 
integral part of these global data sets. 

During the periods 18 August to 4 September 1975 and 2 February to 4 
March 1976, the National Aeronautics and Space Administration (NASA) conducted 
Data Systems Tests (DST-5 and DST-6, respectively) to evaluate several compo- 
nents of the special FGGE observing systems in an operational mode. During 
the DST periods, atmospheric temperature profiles were derived from radio- 
metric data measured by instruments aboard the Nimbus-6 satellite. These in- 
struments were prototypes of those carried aboard the T1R0S-N and NOAA-6 
satellites during the FGGE year. 

As part of the DST evaluations, several research institutions completed 
impact studies designed to assess the effect of Nimbus-6 satellite soundings 
on Northern Hemispheric model predictions. Two and three day forecasts 
started from initial states with and without satellite temperature profile 
data were verified. The results from these studies, as presented by Miyakoda 
et al. (1977) for the Geophysical Fluid Dynamics Laboratory (GFDL) , Hnlen et al. 
(1978) for the Goddard Institute for Space Studies (GISS) and Desmarais et al. 
(1978) for the National Meteorological Center (KMC) , were inconclusive. Con- 
siderable debate was generated between these groups concerning the reasons 
for the lack of significant impact. While the poor quality of satellite sound- 
ings may have been a major cause for the lack of impact, other design features 
of the studies may have contributed. For example, impact was primarily deter- 
mined using two to three day forecasts from relatively coarse mesh (0-400 km 
grid spacing) numerical models. Finer resolution model runs < ver shorter 
periods are known to be more accurate. Also, satellite soundings were mixed 
with data from other sources for these studies. Since satellite soundings 
have different characteristics than radiosonde soundings, care must be taken 
to combine them properly. The results from the DST experiments raised serious 
questions concerning the ability of satellite soundings to provide the • . 
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additlonal, reliable information over data-sparse regions needed for improved 
numerical weather prediction. 

In this paper, a case study approach is used to evaluate Nimbus-6 satel- 
lite soundings over the data-dense region of North America for one day from 
the winter DST-6 period. Analyses of mandatory pressure level heights were . 
constructed on a subset of the Limited-area Fine Mesh (LFM) model grid of NMC, 
which had a grid spacing of 190.5 km. This grid permitted a better resolution 
of synoptic scale patterns than the coarser grids used in the DST impact 
studies. The height fields were constructed using only satellite soundings 
and conventional surface data, eliminating the problems involved with mixing . 
satellite soundings with conventional upper air data. In addition to evalu- 
ating the ability of satellite soundings to define height fields, other inves- 
tigations were undertaken. These included using both the DST resolution and 
a higher horizontal resolution data set. manually checking the soundings for 
horizontal consistency, and subjecting the satellite-derived analyses to the 
LFM model Initialization procedure. 

Several factors have limited the applicability of the results from this 
study to the soundings collected during the FGGE year. This paper presents 
the results from research conducted in 1977 and 1978, a period when the DST 
Impact studies were ^appearing , and before TIROS-N and NOAA-6 soundings became* 
available. The operational retrieval methods used to process TIROS-N sound- 
ings were modified as a result of the experience gained from the DST results. 
Also, the instruments aboard the TIROS-N model satellites have wider scan 
angles than those aboard Nimbus-6, decreasing the wide gaps between adjacent 
satellite passes. Finally and more importantly, measurements in the 15 pra 
infrared channels were degraded by an instrument malfunction during the winter 
DST period, and were not used in the Nimbus-6 sounding retrievals. 

Despite all of these limitations, useful information can be gleaned from 
the results presented. The soundings evaluated in the DST impact tests were de- 
rived in the same manner as those applied in this study. Since satellite sound- 
ings were not mixed with conventional upper air data, certain properties exhib- 
ited by Nimbus-6 soundings that nay have led to the lack of model impact may be 
identified. It is also interesting to compare these "degraded*' results to those 
from more recent studies with TIROS-N soundings such as Schmidt et al. (1981), 
Schlatter (1981) and Phillips ct al. (1979). These comparisons may indicate . 
whether the improvements made in the TIROS-N processing have had a significant 


• effect on the thermal fields defined by satellite soundings. 

The presentation of this study proceeds In the following manner. The 
characteristics of the satellite sounding data are presented in 
Section 2, followed by a description of the study* s design in Section 3. The 
data screening and analysis method appear in Section 4. Section 5 presents 
the results, which are summarized in the conclusions (Section 6). 

2. The satellite sounding data 

The orbital characteristics of the Nimbus-6 satellite, and the retriev- 
al methods used to determine temperatures from the *;easured radiances, played 
an important role in the design of this study. The Nimbus-6 satellite was 
placed in a nearly sun-synchronous polar orbit, with the northbound orbital 
segments passing overhead at local noon (approximately .1800 GMT over the U.S.). 
Consecutive northbound passes cross the equator about 107 minutes apart, with 
passes progressing from east to west. Thus, satellite soundings are available 
in continuous swaths circling the globe, a limited number of which are avail- 
able at the same time as conventional synoptic upper air observations. This 
asynoptic nature of satellite-derived temperature profiles has presented com- 
plex problems in both their application and evaluation. 

Two xnst ruments aboard the Nimbus-6 satel? i ce provided the raw radiance 
data used to derive tropospheric temperature soundings:, the Kigh-resolution 
Infrajred Sounder (HIRS) and the Scan ning Microwave Spectrometer (SCAMS). Both 
instruments had scanning capabilities allowing them to provide fields of data. 
The HIRS scanning geometry Is described in great detail because it served as a 
basis for the method used to convert radiance to temperatures. 

As Nimbus-6 moved along its orbital path, both the HIRS and SCAMS in- 
struments scanned from left to right (see Smith and Woolf, 1976). HIRS had a 
much smaller field of view (29.1 km diameter at nadir) compared to SCAMS 
(144 km). Also, HIRS required frequent inflight calibrations which resulted 
In gaps within the HIRS radiance data called calibration intervals. Each cal- 
ibration interval covered an area equivalent to that of four scan lines and 
was thus approximately 120 km wide in the direction along the orbital path. • 
Twenty scans were made between calibration intervals, each scan consisting 
of 42 individual fields of view. Thus, if a temperature retrieval was made 
for each individual HIRS field of view, 840 soundings would appear between 
calibration intervals, over an area 619 km by 1821 km (roughly equivalent • to* 


one-seventh the area of the adjacent 48 states). The high density of possible 
HIRS soundings is both impractical and unmanageable for most purposes. Also, 
factors such as clouds, varying surface emissivities, pronounced surface ele- 
vation changes, and instrument noise cause some of the radiances to be unsuit- 
able for temperature retrievals. Therefore, HIRS data from several adjoining 
fields were combined to decrease the total number of soundings and to remove 
unsuitable radiances. 

The following description of the method used to combine the radiance 
data in the operational retrieval process is based on the HIRS scanning geom- 
etry illustrated in figure 1. In this schematic diagram individual fields of 
view between two calibration intervals are depicted as circles. (The actual 
field of view would have an elliptical footprint which increases in size as 
the scan angle measured from nadir increases.) 

The first step in obtaining a manageable data set was to subdivide this 
set of 840 fields of view into 30 subset arrays of 28 fields of view each. 

These subsets are called HIRS 4x7 blocks because each contains four fields of 
view along the orbital path by seven fields of view across the orbital path 
(see Fig. 1). Within a 4x7 block, each HIRS measurement is checked for con- 
sistency using SCAMS measurements interpolated to the HIRS fields of view. 
Acceptable radiances are then combined using a method descrioed by Smith and 
Woolf (1976), designed to minimize the effect of cloud contamination. This . 
combination yields a set of clear column (cloud-free) HIRS radiances repre- 
sentative of the 4x7 blocks, located by X's in Fig. 1. These clear column 
radiances are supplemented by SCAMS microwave radiances interpolated to the 
cloud-free radiance locations. 

NESS had originally planned to provide temperature retrievals derived 
at these 4x7 locations during the DST data collection periods. Ar. optimum of 
30 soundings between calibration intervals would then have been available. 
Normal distances between these soundings would have been either 120 km or 
240 km along the orbital path, and 300 kra normal to the orbital path. (The 
longer 240 kn distance along the orbital path is the distance between two 
soundings on either side of a calibration interval.) NMC, however, found the 
number of 4x7 soundings to be too lar^c for operational data handling purposes. 
A maximum of 12 rather than 30 soundings between calibration intervals was 
deemed a more manageable number. To implement this reduction, averaging of 
the clear column infrared and interpolated microwave radiances for adjacent 
4x7 blocks along the. orbital path was performed. 
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Figure 1, A simplified schematic of the HIRS data 
averaging process. Circles represent Individual 
HIRS fields of view. The dashed lines outline HIRS 
4x7 blocks. The 4x7 soun din gs arc indicated by X's 
and the HST soundings by | ] * s . Shaded circles 

represent 4x7 blocks with an inadequate number of 
HIRS measurements for sounding retrievals. 




Figure 2. The LFM grid and the subset 
grid used for the satellite-derived 
analyses. 



Figure^. Nimbus-6 soundings at 
about 'Q600 GMT 22 Feb. 1976. 


i 





-7- 


With only five 4x7 blocks available along the orbital path between cal- 
ibration intervals, overlapping averages were required. In column A (Fig. 1) 
for example, cloud free radiances from blocks 1, 2 and 3 were averaged to 
yield a DST sounding (L]) centered in block 2. Likewise, radiances from 
blocks 3, 4 and 5 were also averaged to create another DST sounding in block 
4. The conversion from averaged radiances to temperature profiles was accom- 
plished with an eigenvector approach also described by Smith and Woolf (1976). 
In an ideal situation, with no unacceptable HIRS data, twelve DST soundings 
would thus have been produced between calibrations with a spacing of 240 km 
along the orbital path within that region, and a spacing of 480 km across the 
calibration interval. The spacing in the other direction would remain at 
300 km. 

. Examples to the right of the orbital path In Fig. 1 illustrate how 
this averaging was affected when adequate acceptable HIRS sounding pairs were 
unavailable within certain 4x7 blocks. These blocks are indicated by shaded 
circles. In column D, two DST soundings result from averaging radiances from 
blocks 2 and 3, and 3, 4 and 5 respectively. Note that the DST sounding is 
positioned at the center of the blocks being averaged. DST soundings also 
appear in blocks El and F3 because only the one block was available for the 
average. . ■ # . 

Both DST and 4x7 temperature soundings are evaluated In this study. 

The roain difference between the 4x7 and DST sounding sets is horizontal reso- 
lution. There are also other differences. In the DST data set, microwave 
only soundings were generated in regions where clear column radiances were 
unavailable, and the microwave data were acceptable. Also, In recalculating 
the 4x7 soundings from the original radiance data, more stringent acceptance 
criteria were used in the clear column infrared radiance determinations. Thus, 
if DST soundings appeared in regions where no 4x7 soundings were made, those 
soundings were either microwave only soundings or ones that passed the less 
stringent DST acceptance requirements. 

3. The study design 

As mentioned earlier a case study approach was used in this investiga- 
tion. DST and 4x7 satellite soundings from five consecutive Ninbus-6 satel- 
lite passes over North America at local noon 22 February 1976, along with 
1800 GMT surface observations, comprise the raw data for the experiment. While 
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the DST data set was cade generally available for the DST impact tests con- 
ducted at several research institutions, the higher resolution 4x7 soundings 
for this February case were generated upon special request by the NESS group 
here at Wisconsin, and are unique to this study. The only previous study eta- 
ploying 4x7 soundings was by Blechman and Horn (1981), who investigated the 
use of higher resolution soundings in delineating a Jet streak over North 
America on 25 August 1975. 

Both the DST and 4x7 data sets were carefully Inspected to detect in- 
consistent soundings, which were removed from these sets yielding two addi- 
tional data sets, the screened DST and screened 4x7 sets. Analyses of heights 
on the mandatory pressure levels were constructed for each of the four result- 
ing sounding sets (both screened and unscreened) on a sv.Dset of the LFM grid. 
These analyses were then passed through the LHl model initialization process. 
Finally, the satellite analysis and initlalizatf -a fields were compared to 
equivalent LFM fields using methods designed *j evaluate both magnitude and 
gradient information. 

Results from these comparisons can provide insight into the following 
pertinent questions. How well do height and temperature fields prepared from 
DST-6 satellite soundings and surface data describe atmospheric features rele- 
vant to numerical prediction? What effect does increasing the horizontal 
resolution of satellite temperature soundings have upon these fields? Are the 
fields considerably improved by a careful manual screening of the temperature 
soundings? And to what degree are satellite thermal analyses affected by .LFM 
model initialization? Answers to these questions, even if they are only par- 
tial answers, could prove valuable in developing better methods for combining 
satellite soundings with data from other sources. The remainder of this r-e';-* 
tion is devoted to explaining the development of this experimental design. 

Common sense dictates that satellite soundings must be evaluated over 
regions with dense conventional data coverage. Since model initializations 
were also tc be evaluated, it seemed logical to employ an existing numerical 
model over a data-dense region for this study. The LFM model fulfilled not 
only these requirements, but also had a 190.5 kn grid spacing (trje at 60°N) 
that allowed better resolution of important synoptic scale features than the 
roughly 400 km grid meshes employed in the DST impact teste. Another advan- 
tage in using the LFM is its familiarity in the meteorological community. The 
analysis subset grid employed in this study is superimposed on the complete 
LFM grid in Figure 2. 
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Several factors Influenced Che choice of cine period studied. A case 
fro. DST-6 vlth strong .inter tenperaturc gradient, would enable a better 
evaluation of satellite capabilities chan a case fro. the summer DST-5. How- 
ever the SIRS Instrument malfunctioned before the start of DST-6. degrading 
e ongwav. Infrared channels end the resultant tenperature profile.. Also, 

cons/T ° r 6 DST ' 6 Peri ° d Wete aVailable here at tbe University cf Wia- 
n in due to a special arrangenont with me. These tape, are not routinely 

rchived and are therefore not readily available for use In case studlee. 

aliy weather maps f r0 . the DST-6 period were studied to find an inter- 
in th! “T? ““ SUlt “ ble “ thiS lnvestl Satlon. A case with two trough. 
The DSt 7 7° " Bl °" appeared fro. 21 February through 23. February 1976. 

from H T tl0 " S eXtr “ Ced £r °" thC NESS - S °“" d1 "* '<*»» 

“ 16 February through 23 February. These locations were plotted over a 

2 °ib’ . l “ 8 " than the Brtd to e “ a 6enc " 1 ld « ° £ '»* <*« 

distribution available. Ninbus-6 passed overhead at local noon and midnight, 

oVcHT "7f ° Ver Bt> “" AnCrlCa “ tC ” '“ >Ura “ elth " *“• ° £ 0M0 GHT and 

ov 2\2 :T aa ° bSa " atl ”" S «" ^ing data readout. 

the , ! “ 0600 GMT - ri ' a “ ltlns ln data -apa over 

the region of Interest, This Is illustrated In Figure 3 for satellite passes 

at roughly 0600 CHT 22 February 1,76. The loch of adegu.te data coverage ' ' 

over the United State, at 0600 CMT precluded any attempt to study the time 

continuity of satellite analyses. It was finally decided to use the „ 

22 F bruary soundings In this case study. The distribution of DST and 6,7 

te sounding locations for this period are shown in Figure 4 

Two properties of the satellite sounding data sets complicated the 

erlf cation procedures used in this evaluation. First, tbe satellite eound- 

mid V ' ' 1 ” aly ' iS rc6 ‘°" “ cco tnketl « roughly 1800 CMT (13 hours), nearly 
dway between the conventional 1200 CMT and 0000 CMT synoptic times, a, 

wn sc ematically in Figure 5. In addition, consecutive passes are 107 

« nutes .part, which leads to a 6 hour time difference between the eastern- 

“ 7 " aS ; an, " 0St S< “' ,nd1 " 88 £ " '*» region. This problem 1. better ' 

da a .7' * r EUr0 6 ' from the unscreened DST 

l ea ' T “ P " blCm 1S ” 1S ° da,,,CtCd thlS ««•"• "a shading de- 
that a 7 Sat0, “ te regions between consecutive passes 

ppt.ar sout of about 55°N. Those data gi . ps cause unique analysis 
problems. 
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Figure 4. ■ Ninbus-6 satellite soundings around 1800 GMT 
22 Feb. 157.6,' Panel A - DST sounding locations, panel B - 
4x7 sounding locations and panel C - final 4x7 sounding 
locations .supplemented with qualifying DST soundings (see 
text). □ indicates soundings removed in the screening 
process . 
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SATELLITE SOUNDING OBSERVATION TIMES 

Figure 3. A schematic illustrating the time interval of the 
satellite soundings in relation to the four synoptic, times 
used to derive the verification data. 



Figure b. An analysis ot satellite sounding observation times. 
The ^aps between consecutive satellite passes are shaded. 
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Rather than just comparing the satellite results at roughly 1800 GMT 
to the bracketing 1200 GMT and 0000 GMT values, as done by Horn et al. (1976) 
and Petersen and Horn (1977), bogus LFM verification fields were generated 
which incorporated these limiting features of the satellite sounding data sets, 
namely data gaps and asynoptic times. The first step in constructing bogus 
LFM analyses was to remove the time differences between the satellite and LFM 
analyses. This was performed by interpolation to satellite sounding times 
from the four bracketing LFM analyses using overlapping quadratic polynomials. 
Two sets of bogus LFM analyses were constructed. The first simply interpo- 
lated grid point values in time using the grid point analysis of satellite 
observation time illustrated in Figure 6. The generation of the second bogus 
set involved a more complicated two step procedure. In the first step, grid 
point values were interpolated in space to the satellite observation locations 
for all four bracketing LFM analyses, and then interpolated to the reported 
satellite sounding time. (These soundings were used as the verification data 
in the colocated comparisons.) The bogus LFM set was then derived by inter-' 
pointing these colocated values in space back to the grid points using the 
same analysis and filtering procedure applied in the satellite sounding thick- 
ness analyses. 

The first bogus LFM analysis set derived from time interpolation only 
provides an estimate of the LFM fields valid at the times of the satellite 
observations. In the second bogus LFM set, both time and space interpolations 
are employed to further simulate inherent spatial limitations of the satellite ' 
soundings in this case. This latter bogus LFM set is labelled INTERP LFM to 
emphasize the spatial interpolation. A comparison between these two bogus LFM 
analysis sets given later in this report provides an estimate of how much in- 
formation is lost in the analysis process. 

The satellite-derived analyses can at best be expected to reproduce 
the INTERP LFM bogus analyses. Therefore, those fields serve as the standard 
of comparison in the visual and gridded statistical evaluations. In essence, 
the LFM analyses have been degraded to the same level as their satellite- * • 
derived counterparts to provide equivalent comparisons. 

A final point to consider in this discussion of the case study design 
is that the Nimbus-6 soundings provided only temperature data at 21 levels 
between 1000 mb and 100 mb inclusive, and thickness values derived from these • 
temperatures for layers between 1000 mb and each remaining mandatory level 
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(1000-850 mb, 1000-700 mb, etc.)* Reference level Information was needed to 
determine heights on mandatory pressure levels used in the evaluations. 

Heights of the 1000 mb level derived from surface observations taken at 
1800 GMT from the surface synoptic stations and other sources provided the 
reference level information. These 1000 mb heights were derived using stan- 
dard reduction methods described in SchUepp et al. (1964). 

In constructing the bogus LFM data sets, thickness values relative to 
1000 mb rather than heights were the fields interpolated in time. The same 
reference level values used in the satellite analysis were used to derive the 
heights in these bogus LFM sets. The reason for proceeding in this manner, 
rather than interpolating height values directly, was to remove any bias that 
may enter from the different 1000 mb height determination methods used in this 
study and at NMC. 


4 • Data screening and analysis methods 

An important and time consuming segment of this study was the manual 
screening of the DST and 4x7 satellite soundings. This screening procedure 
involved a careful examination of both vertical and horizontal sections of 
the atmosphere to detect inconsistent soundings removed to create the screened 
data sets. Before proceeding with this screening, the original 4x7 data set 
was supplemented with selected DST soundings in data gaps such as Texas, the 
northeast U.S. and central Canada, where soundings available in the DST set 
were not available in the 4x7 set (see Figures 4A and 4B) . Only 14 DST sound- 
ings derived from a single 4x7 block of HIRS measurements were added to the 
original set of 356 4x7 soundings. Those additional soundings were either 
mi wave only or those rejected due to more stringent 4x7 acceptance criteria. 
While such soundings may have been of lower quality than the original 4x7 set, 
they provided more information than no data in these regions, and still had 
to pass the consistency checks in the screening procedure. The positions of 
the amended 4x7 sounding set, hereafter called simply the 4x7 set, are shown 
in Figure 4C. 

The first step in the screening procedure was to plot and hand analyze 
horizontal maps of mandatory level temperatures and thicknesses relative to 
1000 mb from both the DST and 4x7 sounding sets. Inconsistent soundings were • • 
noted for each analysis. Cross-sectional analyses of potential temperature 
and isotachs of the thermal wind relative to 1000 mb were constructed from 
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the level ter.pcratu^e information using an objective method developed by 
Whittaker and Petersen (1977). Thermal wind analyses measure thermal gradi- 
ents along the cross section and are valuable in detecting vertically inte- 
grated gradient errors. Six roughly north-south cross sections were con- 
structed for each orbital path, one for each of the columns labelled A through 
F in Fig. 1. The 6 cross sections for each pass were displayed sequentially 
to aid in detecting inconsistencies from one section to the next. . : 

Figure 7 illustrates how cross-sectional analysis can emphasize sound- 
ing errors. The westernmost column of soundings from the 4x7 data set were 
used in this figure. Often errors of one sign in the lower troposphere for 
a given sounding were accompanied by errors of the opposite sign in the upper 
troposphere. The two soundings labelled 5068 and 5090 show this typical probr 
lem, with 5068 being too warm in the lower troposphere and too cold near 
300 mb, while the opposite is true of 5090. Note how the thermal wind analy- 
sis emphasizes a problem between 5084 and 5090. 

Inconsistent soundings detected in the horizontal hand analyses were 
double-checked on the cross sections, with gross inconsistencies at any level' 
causing removal of the entire sounding. The multidimensional view from both 
horizontal ar.d vertical sections helped in making the decision of which were 
the poor soundings, , a decision that was sometimes easy, but often difficult.. . 
In the example from Figure 7, the decision to remove 5068 and 5090 from the 
screened data set was fairly straightforward, but the decision to remove 5084. 
was more difficult because it required comparisons with other cross sections. 

The soundings deemed inconsistent by this method are indicated by a 
□ In Figure 4. The yield of acceptable soundings is much higher for the 
DST set than for the 4x7 set. Only 16 of 300 (5.32) of the DST soundings 
were removed to form the screened DST set, while 107 out of 370 (28. 9X) were 
removed from the 4x7 set. This result is not surprising since as many as 
three 4x7 soundings were averaged together to form one DST sounding, smooth- 
ing errors found in the 4x7 set. Thus, while the yield of the DST soundings 
was higher than that of the 4x7 set, good data may have been averaged with 
poor data to give mediocre yet consistent soundings. 

Unfortunately, cloud contamination is not the sole factor in producing 
inconsistent soundings. A distinct inconsistent-cloudy relationship is dif- 
ficult to perceive from Figure 8, where the 4x7 sounding locations are super- 
imposed upon a hand analysis of percent cloudiness reported for each satellite 
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Figure 7. A sample cross section used in the screening procedure* 
Isentropes (°K) are solid, and isotachs of the thermal wind com- 
ponent relative to 1000 mb and normal to the cross section (m s“l) 
are dashe.d. Soundings 5090 and 5084 are 120 km apart. 
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Figure 8. Satellite sounding percent cloudiness analysis 
superimposed upon the 4x7 sounding locations (Fig. 4C) . 
Isopleths are draw every 25%, with areas above 50% cloudi- 
ness shaded. [7^3 again indicates unacceptable soundings. 



- 10 - 


sounding. (Comparisons against a cloud analysis from the surface synoptic 
network, satellite cloud photographs, radar maps and precipitation reports 
demonstrated that this analysis gives a reasonable representation of cloud 
cover,) The cloud band along the east coast was by far the most active in 
terms of precipitation, which may have degraded microwave only soundings in 
that region. Much of the cloud contamination problem already appeared in the 
form of missing data in the mostly cloudy regions. 

The satellite data screening procedure took considerable time ond ef- 
fort. The decision to remove a given sounding often required a complete, 
three-dimensional view of the data structure. This fact, along with the lack 
of any systematic relation between clouds (and precipitation) and inconsistent 
soundings, would complicate attempts to computerize this screening procedure. 

The next step in the study was to construct gridded height analyses for 
both the screened and unscreened 4x7 and DST sounding sets (4 satellite- 
derived sets total). The process of transforming data at observations into 
final analyses on the LFM grid involves several interpolation and filtering 
operations which form the analysis procedure. Several characteristics of the 
basic data (satellite soundings and surface observations) presented certain 
difficulties which required special consideration in the development of this 
procedure. 

The wedge-shaped data gaps between satellite pass south of about 55 N 
create an abrupt transition in data coverage that can cause severe interpola- 
tion problems. Also, the satellite soundings provided only thickness infor- 
mation which was combined with base level information from surface observa 
tions to form the isobaric height data used in the analyses. The analysis 
procedure was designed taking these factors into account* 

The analysis process employs two types of interpolation: interpola- 

tion from the unevenly spaced observations to the uniform set of grid points , 
and its inverse, interpolation from grid points to observations. The latter 
is needed to provide first guess and 1000 mb height values nt observation 
locations. Several methods were tested for each tvpe of interpolation, using 
a typical set of Nimbus-6 sounding locations to simulate the data gap problem, 
Data values at both I.VM grid points and those observation locations were spe- 
cified by analytic functions designed to approximate atmospheric height ind 
temperature fields. The ability to specify values at both observation and 
grid point locations permitted an exact measure of accuracy for the different 
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interpolation methods tested. 

The results from these interpolation experiments arc presented in 
Koehler (1979). The grid to observation interpolation test results indicated 
that the optimum interpolation method developed by Gandin (1963) was among 
the most effective for this particular application. A modification of this 
method also discussed by Gandin, called optimum interpolation with normalized 
weights, was chosen for the observation to grid interpolations in this study. 
Normalization of the interpolation weights was incorporated to counter the 
variability in the sum of the interpolation weights that appeared in the 
satellite data gaps* 

Excellent discussion on the theoretical basis for optimum interpola- ^ 
tion methods and their development are presented by Gandin (1963), Schlatter*. 
(1975) and Bergman (1978) . A detailed description of the interpolation method 
used in this study is presented in Koehler (1979). Only certain features of 
this method will be discussed he.e. 

The method. is designed to incorporate the statistical structure of the 
field being analyzed (isobaric height). The observational error structure is 
also included in the analysis method, such as the correlated nature of the 
satellite soundings. The specification of these statistical properties is 
fully described in Kdehler (1979). 

The eight closest observations are used for a given grid point, and 
only one set of interpolation weights are determined for each data set. These 
weights are used to construct analyses at all mandatory lcvelh. The analyses 
for most variables used persistence first guess fields taken from the LFM 
analysis 18 hours previous to the satellite analyses. Analyses of variables 
with no first guess field (such as the time analysis in Figure 6) could still 
be performed due to normalization of the interpolation weights. 

As mentioned before, a grid to observation interpolation method was 
needed to transfer the gridded first guess values and 1000 mb heights to the 
satellite observation locations. A method introduced by Bleck and Haagenson 
(1968) that employs overlapping quadratic polynomials proved to be both ac- 
curate and efficient for this type of interpolation. 

A step-by-step description of the analysis procedure follows. Since 
all four satellite sounding sets and the bogus LFM fields used the same 
1000 mb base level information, the 1800 GMT 1000 mb analysis was completed 
fj.rst. First guess values were interpolated to the surface synoptic station 
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antf ship report locations. Interpolation weights were then derived for the 
8 closest observations to each grid point, with only those observations within 
2000 km eligible for use at that point. Tnus, it is possible that fewer than 
, . 8 observations were used in the interpolation to certain grid points in data 

sparse regions. After the analysis values were determined by applying these 
weights, a filter consisting of a smoother-desraoother was passed over the 
\ 1000 mb analysis to remove small scale noise from the initial analysis. This 

filter is the same as one applied in the postprocessing of LFM initial hour 
and forecast fields. (See Gerrity and Newell, 1976.) The spectral response 
of this filter, and another called the combined filter are depicted in 
Figure 9, the latter being used to filter the upper level, satellite-derived 
height analyses. Note that the effect of these filters is to remove all 2Ax 
noise from the analysis, while synoptic scale features at longer wavelengths 
are retained. 

The next series of steps was repeated for each of the four satellite 
sounding sets. The 1000 mb height values at each satellite sounding location 
were obtained by* a grid to observation interpolation from the 1000 mb analy- 
sis. These 1000 mb heights were added to the satellite-derived thicknesses, 
yielding isobaric heights at each of the nine mandatory levels between 850 
and 100 mb. First guess values were also estimated at . each sounding location 
by grid to observation interpolation. 

Analysis weights were then calculated for each of the four sounding 
\ location distributions. The raw height analyses at all mandatory levels for 
a given sounding distribution were then calculated with one set of interpola- 
tion weights. The combined smoother from Gerrity and Newell (1976) was ap- 
plied to these raw height analyses to produce the final satellite height 
analyses studied extensively in the following sections. 

One facet of the experiment was to determine the effect of model ini- 
tialization on the mass field variables described by satellite soundings. The 
LFM model initialization program at NMC, described by Gerrity (1977), was run 
with the satellite data sets. While this initialization consists of several 
operations, the vertical interpolation of height and temperature values from 
the 10 mandatory pressure levels between 1000 mb and 100 mb to the model s 
o-coordinate layers is of primary importance here. Also, while values cn o 
surfaces serve as model init ializa t ion fields, researchers using LFM save 
tapes are normally restricted to studying postproeessed initial hour data on 
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Figure 9, Filter response versus wavelength, (See Gerrity 
and Newell, 1976.) 


Table 1. Statistics for the second bogus LFM minus first bogus LFM 
height field differences. Bias and RMS values are in meters while 
SI values are dimensionless. These statistics provide a measure of 
the information lost in the analysis process due mainly to systematic 
data gaps between satellite passes. 
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isobaric surfaces. Another vertical interpolation from o to p is involved 
in the postprocessing, which has the effect of further smoothing the height 
and temperature field (Hyde, 1977). This should be kept in mind when con- 
sidering the effects of the initialization presented in subsequent sections 


5 , Results 

The degree to which satellite soundings can provide either magnituda 
or gradient information for atmospheric mass field variables is an important 
consideration in their future application. Several methods were employed to 
evaluate both magnitudes and gradients from the satellite-derived analyses.. 
These methods include visual comparisons of the basic fields, statistical . 
measures commonly used in forecast verifications and comparisons of dynamic', 
quantities derived from the basic fields, such as geostrophic wind, geo- 
strophic temperature advection and available potential energy. 

As mentioned earlier two LFM verification data sets were constructed 
from a time interpolation of the bracketing LFM fields to the times of the. 
satellite observations. One of these sets (the 1NTERP LFM) included an ap- 
plication of the horizontal analysis procedure, while the other did not. This 
section begins with a comparison between these two bogus LFM data sets, which 
can be used to estimate the effect that data gaps between satellite passes 
have upon the satellite analyses. This comparison is followed by the evalu- 
ations of the satellite soundings described above. 

A complete and detailed discussion of the results using data from 4 
satellite sounding sets at 10 mandatory levels for both analysis and ini- 
tialization fields could become quite laborious. Many results are given in 
tabular form for all mandatory levels. However, the text devoted to these 
tables will usually be brief and will concentrate on points of greatest 
relevance and interest- 

In many of these presentations, labels for the 4 satellite sound g 
sets have been abbreviated. The unscreened DST and 4x7 sets arc termed ALL 
DST and ALL 4x7 respectively. Similarly, their screened counterparts arc 
labelled SCR DST and SCR 4x7. 

a. A comparison between the two bogus LFM data sets 

The first, bogus LFM analysis set derived only from time interpolation 
provides an estimate of the LFM fields valid at the times of the satellite 
observations. This is illustrated in Figure 10 which includes heights. 
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temperatures, and isotachs of the geostrophic wind ar J00 mb for the first 
bogus LFM and bracketing 1200 GMT and 0000 GMT LFM data sets. (Diagrams in 
this and subsequent figures valid at satellite observation tit.es, are labeled 
1800 GMT.) As indicated in Figure 10, the general propagation of features 
such as troughs, ridges and jet streak patterns, are well represented by the 
time interpolation. A slight deg -adution of the height gradients can be de- 
tected in a comparison of the maximum speeds for the jet streak in the south- 
east United States. The speed decreased from the bracketing values of 85.8 
and 89.8 m s' 1 to a value of 83.5 m s -1 . 

In the second bogus LFM set, both time and space interpolations were 
employed to further simulate inherent spatial liraitati.ns of the satellite 
soundings. A comparison between the first and second bogus LFM analyses is 
offered in Figure 11. The second bogus LFM set is labeled INTERP LFM to em- 
phasize the space interpolation. Gaps between satellite passes are shaded. 

The largest height differences of over 40 m are found in these data 
gaps. Unfortunately, part of the jet streak in the southeast is also located 
in one of the data gaps and suffers a further decrease to 71.1 m a -1 , due to 
space interpolation and smoothing of the height field. Bias and RMS differ- 
ences, and SI scores for the second bogus LFM set versus the first are pre- 
sented in Table 1. While the bias and RMS differences increase with eleva- 
tion up to the tropopause, the SI score which is designed to measure gradient 
difference decreases. The SI score of 14.8 at 300 mb corresponds well with 
the 14.92 decrease in the jet max value, which suggests that gradient losses 
are distributed over the entire analysis region. Such a comparison provides 
an important measure of the loss of information due to systematic data gaps 
between satellite passes and the analysis procedure. Remember, satellite 
analyses can at best be expected to reproduce this second LFM bogus set. 

These fields will then serve as the standard of comparison for the visual and 
** t icdl evaluation methods. 

b. Layer Dean temperature comparisons at observation locations 

Following an approach used by Hayden (1977), mean temperatures for' 
layers between consecutive mandatory pressures were calculated at observation 
sites, and compared to similar values interpolated from gridded analyses to 
the observation locations. Comparisons with the four satel'itc sounding sets 
were made against both bogus LFM and satellite-derived analysis values inter- 
polated back to the sounding locations. The former comparisons against LFM 
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Figure 11. Bogus LFM 300 mb comparisons. Panel A - height (solid, dam) 
and temperature (dashed, r C) analyses for the first bogus LFM set (time 
interpolation at ; rid points only). Panel C - same as A except for the 
final bogus LFM set (time interpolation at sounding locations and sub- 
sequent space interpolation to grid points). Panels B and D are the 
corresponding geostrophic isotach analyses (in s“l). Panel E shows the’ 
height difference fields (dan) between A and C. 
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analysis values measure the accuracy of the soundings, while the latter com- 
parisons against satellite-derived analyses measure the noise removed from 
the original data by the analysis process. Included in the discussion are 
Hayden’s comparisons of radiosonde and Nimbus-6 soundings against NMC'a 
Northern Hemispheric analyses for the period from 15 February to 20 February 
during DST-6. While Hayden’s comparisons were hemispheric, comparisons for 
the DST and 4x7 data sets used in this study were restricted to sounding lo- 
cations within the satellite analysis subset of the LFM grid. Bias and stan- 
dard deviation statistics for the observation minus analysis temperature dif- 
ferences are given in Table 2. 

Values from Hayden's RAOB vs. analysis and Nimbus-6 vs. analysis com- 
parisons arc shown first. The hemispheric analyses used in his comparisons 
are based primarily on radiosonde measurements and a^e independent of Nimbus-6 
satellite observations. This explains tb~ relatively small bias differences 
in the RAOB vs. analysis comparisons in Table 2. The standard deviation re- 
sults were larger for the Nimbus-6 satellite soundings than for the radio- 
sonde measurements, especially near the surface and tropopause. Bia 9 dif- 
ferences in the Nimbus-6 set arc smaller below 700 mb than those from the 
radiosonde set, which may indicate a problem in hemispheric analyses at lower 
elevations. 

The ALL DST set from this study Is equivalent to Hayden's Nimbus-6 
sounding set in terms of the sounding retrieval and averaging techniques de- 
scribed earlier. Standard deviations for the ALL DST vs. LFM analysis dif- # 
ferenccs are larger than those from the Nirabus-6 comparison, again at lower 
levels and the tropopause. The ALL DST bias differences arc also larger than 
the Nirabus-6 biases at lower levels, and the maximum bias of 2.44"C for the 
ALL DST set is in the 250-200 mb layer. Several factors contribute to dif- 
ferences in these ALL DST and Nimbus-6 comparisons (Table 2). The ALL DST 
results arc only for a limited portion of the LFM grid over North America 
on February 22, while Hayden's Nlmbus-6 comparison was hemispheric for five 
days during an earlier period. Also, the verifying analyses for the two com- 
parisons (LFM vs. hemispheric) were constructed with different analysis 
methods on grids of different resolution. 

The major difference in the comparisons of the DST and 4x7 sounding 
sets against analysis values in Table 2, is that bias differences for the DST 
sets are larger near the sutface and smaller near the tropopause than bia-.sB. 
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Table 2. Layer mean temperature comparisons (°C) between observed 
values and analysis values Interpolated to the oo:*ervation loca- 
tions. The RAOR vs. Analysis and Nimbus-6 vs. Analysis comparisons 
are from Hayden (1977). 


BIAS 


IAYIR (MB) 


. J 

1000 

TO 

B30 

nso 

TO 

700 

700 

TO 

300 

500 

TO 

400 

4 00 
TO 
500 

300 
TO 
2 30 

230 

TO 

200 

200 

TO 

130 

130 

TO 

too 


OAT A 


SI IS 


NAOB VS • ANAL 
NIMBUS - 6 VS ANAL 

on 
0 3 

•02 
-0 1 

0 2 
0 3 

0 0 
03 

OO 

03 

06 

1 4 

*0 2 
06 

03 

06 

09 

03 

AIL 

OST 

VS 

If M 

anal 

1 II 

1 60 

096 

1 21 

051 

1 06 

2 44 

0 71 

0 09 

SCR 

DS* 

VS 

LtM 

ANAL 

120 

1 62 

09b 

1 21 

033 

1 09 

242 

073 

006 

ALL 

4X7 

vs 

If M 

ANAL 

0 74 

1 10 

0 20 

0 33 

003 | 

| 1 36 

303 

102 

0 02 

SCH 

4X7 

vs 

l 1 M 

ANAL 

1 01 

1 20 

0 2ft 

031 

0 00 j 

[ 1 24 | 

2 66 

0 03 

006 

Al L 

osi 

VS 

•SAT 

AN Al 

0 01 

0 09 

0 05 

-001 

-003 j 

I 0 01 

001 

0 02 

0 02 

SCR 

OST 

vs 

SAT 

ANAL 

0 05 

OOf 

004 

-0 01 

j *0 02 | 

i 0 02 

O 01 

000 

0 00 

Al l 

4x r 

vs 

S Al 

ANAL 

0 20 

0 35 

0 25 

0 12 

1 ~ j 

i 0 04 

i 

- 0 16 

*0 22 

j * O 1 3 

*002 

SCH 

AX 7 

vs 

SAT 

AN Al 

O 14 

0 12 

0 0 « 

O 02 

0 06 

-002 : 

-O 06 

j -002 

001 


SFANOAKO OfVlATlON 






t AUN 

IM0) 




1000 

630 

700 1 300 

400 | 

g 

r, 

• j 
6 

130 


— . . * 

TO 

TO 

TO j TO 

TO | 

TO | TO t TO 

TO 


PA1A SLTS 

030 

700 

■>00 i 400 

300 

230 I 200 1 130 

too 

; 

RAUO VS ANAL 

1 6 

TV- 

13 ! 13 

,f i 

2 0 { 17 ! lb 

1 9 


NIMBUS- 6 VS AN Al 

2 b 

2 0 

♦ 

1 7 j 17 

•» 

23 ; 31 ! 23 

2 0 


All OSf VS l » M ANAL 

3 nr 

2 4 7 i 16 3; 1 b9 

1 bb 

3 10 4 19 , 2 49 

f ‘ - 

1 73 


SCR OST VS It M AN Al 

X R9 

2 99 

101 ; 167 

1 bb 

3 13 ! 4 26 j 2 30 

i . 

1 72 


ALL 4X 7 VS l t M AN Al 

\ 

3 ;o 

2 6 7 I 1 64 l 173 

1 64 | 2 94 i 401 | 231 

1 70 


SCH 4X 7 VS ItM ANAL 

3 60 

t * 

2 6b | 106] 1 69 

1 67 

2 40 1 396 | 2 47 

1 39 


AIL OSJ VS SAT ANAL 

1 4 3 

O *4 

0 7b j 0 63 

0 72 

0 74 | 0 73 j 0 63 

0 70 

/ 

SCH OST VS SAT AN Al 

1 39 

ObB 

0 bb 1 0 68 

0 37 

0 b 3 , 0 69 j 062 

073 

/ 

All 4X7 VS SAT ANAL 

t 42 

t 90 

1 39 j » 34 

1 0 7 

l 43 j * 36 ; 118 

1 21 


SCR 4X7 VS SAT ANAL 

1 32 

|iwr 

0 74 | 0 73 

i 0 bb 

0 74 ! O 6 4 | O 63 

0 66 


Nl'MtUH 01 SOeNOtNuS 
All OST * 104 All 4*7 - 222 

SCR OST ' 17** MP 4X7 - 133 








■w w«*i r. W^rof<f!W^yTOT> J>V->,K'W V ' J ."V » v«hj 

-26- 






from the 4x7 sets. Little or no improvement from screening the satellite data 
sets can be detected in these comparisons against LFM analysis values. 

Comparisons of satellite observations versus satellite-derived analysis 
values provide more information. While a general label of SAT ANAL is used in 
Table 2, a given satellite sounding set was only compared to analyses con- 
structed from that set. For example, the screened 4x7 soundings (SCR 4x7) 
were only compared with SCR 4x7 analysis values. The small bias differences 
indicate that the analysis procedure does preserve the means of the data 
fields. Standard deviation values for the ALL 4x7 case are almost twice those 
for the ALL DST case, indicating more frequent inconsistencies in the raw 4x7 
data set. The screening procedure succeeds in removing inconsistent soundings 
as shown by smaller standard deviations in the screened cases. The effect of 
screening the ALL 4x7 set is substantial, although the standard deviations for 
the SCR 4x7 vs. SAT ANAL comparisons are still slightly higher than those from 
the SCR DST vs. SAT ANAL comparisons. 

c. Height and temperature analysis comparisons 

The effect of a high noise level in the raw data on resulting analyses 
is illustrated in Figure 12 for the ALL 4x7 sounding set. Anomalous troughing 
in Colorado, Florida and the extreme west central portion of the grid stands 
out in the 500 mb analysis. All three problem ureas are in data gaps and re- 
sult from noisy soundings. The detrimental effect of a noisy sounding near 
the center of a satellite pass will be moderated by surrounding observations 
during the analysis. However, the analysis method tends to extrapolate gra- 
dient information into data gaps. Any noise near the edges of a satellite 
pass will create anomalous gradients which are extended into the data gaps. 

For example, heights for the ALL 4x7 sounding in central Colorado (Figure 4C) 
were too low, which caused an anomalous gradient between this sounding and 
the one Just to the east on the Kansas-Nehrnskn border. The extrapolation 
of this gradient into the data gap formed the anomalous trough in Colorado. 
Because of these unsatisfactory results the ALL 4x7 data was not used in pre- 
paring model initializations or in performing more statistical evaluations. 

Height and temperature analyses for the three remaining satellite data 
sets are presented In Figures 13 through 16 at 850 mb, 500 mb, 300 tnb and 
200 mb respectively. The final bogus LFM analyses are included for compari- 
son. Trovgh and ridge positions are generally well represented in the 









satellite-derived analyses, however, as might be expected, the three satellite 
sets resemble each other more than they resemble the bogus LFM data. This is 
not surprising since all three satellite sounding sets were derived from the 
same HIRS and SCAMS measurements. In fact, the subtle differences between the 
ALL DST and SCR DST analyses are almost imperceptible in these figures. 

At 850 mb (Figure 13) , thermal gradients in the northeast corner of the 
grid are much weaker for the satellite cases than for the LFM, This is due to 
the inability of the satellite measurements to capture the colder temperatures 
north of the Great Lakes, Assuming that the LFM's temperature values over 
that portion of Canada to be correct, the satellite-derived temperatures are • 
almost 7°C too warm. (The strength of this 850 mb cold dome in the LFM analy- 
sis is supported in the Canadian radiosonde reports in that region at both 

» 

1200 GMT February 22 and 0000 GMT February 23.) Note also that the warm re-= * 
gion indicated by the 4°C line in Wyoming and Montana in the bogus LFM analy- 
sis is not well defined in the satellite-derived analyses, Indicating a slight 
cold bias in that region. 

The temperature field at 500 mb (Figure 14) is well represented in. the 
satellite analyses. However, the overestimation of 850 mb temperatures Is re- 
flected by higher satellite heights in the 500 mb trough. A short wave trough 
not found in the LFM analysis is indicated in the satellite data in associa- 
tion with an Alberta low present in the conventional surface analysis (not 
shown). It is possible that the widely spaced Canadian radiosonde network 
was unable to define this feature. If so, the ability of satellite soundings 
to better define such shortwave features over Canada would be a definite ad- 
vantage. On the other hand, in another data-sparse area, the trough south of 
California in the satellite analyses is not supported in the available conven- 
tional data. 

The inability of satellite soundings to define the 300 mb cold pool 
over the Great Plains (see Figure 15) is an indication of their known diffi- 
culties in defining the temperature structure near the tropopause. This Is 
further accentuated at 200 mb (Figure 16) where satellite soundings were able 
to locate a cold pool over the western United States, as indicated by the 
-56°C line. However, the intensity of this cold pool is underestimated, and 
subsequently, the strong reversal of the temperature gradient evident at this 
level in the LFM data is greatly underestimated. 

A set of height difference fields (satellite minus LFM) at 300 mb is 
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Figure 17. Height difference fields (dam) between the satellite-derived analyi 
LFM analysis at 300 mb. Negative differences arc dashed. 
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shovn in Figure 17. Positive values in this figure indicate areas where 
satellite-derived heights were higher than the LFM heights. The general pat- 
terns for the three satellite sounding sets are similar, with the greatest 
overestloatlons In troughs and smaller underestimations in ridges. Positive- 
120 a differences In the eastern trough correspond to 1000 mb to 300 mb layer 
mean temperature differences being 3.4*C too warm, due mainly to low level 
temperature overestimates such as the 7"C error at 850 mb mentioned earlier. 
An anomalous trough in the SCR 4x7 set is evident fro* the large negative 
differences south of Florida, which will affect several of the statistical 
measures of accuracy. 

The lack of conventional data in the Pacific region of the analyses •* 
makes verification there less reliable. To focus more attention on the land 
areas, the analysis region was divided into Pacific, western North American 
and eastern North American subregions, as shown In Figure 18. This also 
separates the flow pattern regimes, with the Eastern and Pacific regions 
dominated by troughs, and the Western dominated by a ridge. .. * 



Figure IS. Subdivision of the analusia region for the calculation of 
evaluation ciat ictics. 

Height statistics for the three subregions, and the entire region, are 
given in Table 3. Positive bias differences predominate except in the tropo- 
sphere of the Western region, where small negative values appear in the ridge 
regime. Both bias and RMS differences increase with altitude in all regions, 
while SI scores decrease from 700 mb up to 300 mb in all except the Pacific 
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region. A closer examination of how the SI score is fontjlated will help ex- 
plain the decrease in this score with altitude, up to the 300 mb level. The 
SI score is defined as the sua of the magnitudes of satellite gradient minus 
verification (LFM) gradient vector differences, divided by the sum of the ver- 
ification (LFM) gradient magnitudes. As altitude increases in the troposphere, 
LFM height gradients generally increase in magnitude, while gradient differ- 
ences between the satellite and LFM fields also increase, but at a slower rate 
than the LFM height gradients. In terms of the SI formulation, the denomina- 
tor is increasing faster than the numerator, yielding decreasing SI scores. 

The anomalous trough in the southeast corner of the SCR 4x7 Analyses 
produces smaller biases than the DST sets in the Eastern region. However, 
this erroneous trough is characterized by larger RMS and SI scores which ac- 
company the smaller biases. In fact, the SCR 4x7 heights exhibit higher SI 
scores everywhere except at 700 and 500 mb in the Western region. In both 
the Western and Eastern subregions, smaller biases accompany larger RMS dif- 
ferences for the SCR 4x7 aet compared to the DST sets. Only in the Pacific 
do both smaller SCR 4x7 bias and RMS differences appear together. Unfortu- 
nately, verification is questionable in that region. 

Layer mean temperature comparisons are presented n Table 4. Through- 
out the troposphere, satellite estimates are too warn, as indicated by the 
positive biases-. Diases of one sign at lower levels are not being compen- 
sated by biases of the opposite sign at higher levels, as found in studies 
with other satellite data sets such as Horn et al. (1976). If such compen- 
sating biases exist, an advantage can be found in using height and thickness 
variablca in which temperature Information is integrated vertically. This 
allows biases of opposite sign to counterbalance each other. This compensa- 
tion is not present in DST-6 soundings, and thus, the advantage in using 
vertically-integrated variables is lost. 

While RMS values f or raid and upper tropospheric layers (700 mb tc 
300 mb) are reasonable, values near the surface and tropopause are unreanon- 
ably large. The strong positive biases and related large RMS values at lower 
levels in the Eastern trough region have already been alluded to in the 850 mb 
analysis, and 300 nb height difference field discussions. Very large SI scores 
of over 100 are found in the tropopause layers and are indicative of the tem- 
perature gradient reversal problem also described earlier. 

Generally , the satellite data’s deficiencies in defining the themal 



same format as Table 3. In this case 


WESTER*! REGION 








wm 

H 


lisa 

151 

!B51 



rrn 


ta 

H 

s 

m 

H 

H 

B 

gaga 


IEQ 



eg 

KO 


m 

CD 



EES 


iS 


IBS 

ipw 

m 

m 



Bg 


DB 

ICT 

fff 

Bp 

m 

B D 


BB 

BB 


BO 

fffl 

na 


m 

ra 



ED 


ED 

B3 

m 

QQ 


m 


na 

Bn 


B3 

KEB 

m 

m 


mg 


ini 

m 


m 

m 

m 

rm 

lH 

ED 

■EBBS 

Bill 

m 


na 

ED 

m 

m 


m 


BUI 

BUI 


m 

ESI 

ED 


CD 

BD 


EASTERN REGION 




































is- 


structure in the cold dome account to a large degree for problems in deflnins 
the tropospheric height fields. As expected, bias and RMS mean tempers Hue 
values for the whole region analysis comparisons arc similar to those from 
the raw data comparisons at observation points presented earlier. Only slight 
improvements can be detected in the RMS differences due to screening the DST 
data set. 


d. Energetic and dynamic parameters * 

Parameters related to atmospheric energetics and dynamics derived di- 
rectly from mass field variables (height, temperature and pressure) can as- 
sist in the satellite data evaluations. The first parameter to be discussed 
is available potential energy (APE), a quantity that provides a one number 
measure of the integrated bnroclinlty over a region, and has been used in 
other satellite data evaluations (Pesmarals et al., 1978). A technique de- 
scribed by Koehler (1979) was employed in the APE calculations. It uses an 
exact formulation of APE in lscntroplc coordinates rather than an approximate 
formulation in tsobaric coordinates used in previous studies. 

The available potential energy values presented in Table 5 are in tha 
form of specific (per unit mass) values. The mass in the computation voluma 
can vary for different data sets, so effects of the' total mass in the APE val- 
ues have been removed. The calculations were made with an upper boundary at 

the 370K lscntroplc level, which corresponds roughly to the 150 mb isobarlc 
surface. 
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The strong baroclinity associated with the east coast trough is re- 
flected by the large APE values for the Eastern subregion. Satellite esti- 
mates show about a 10Z Weaker APE in calculations for the Eastern region com- 
pared to LFM values. Over the entire region, APE is 18% weaker for both DST 
analysis sets than for the LFM set. APE values calculated from the SCR 4x7 
analyses are slightly larger than those from the DST analyses, and are thus 
closer to the LFM value . This normally implies higher baroclinity in that 
set than in the DST data sets. Without other gradient measures such as the 
SI score, it would be difficult to say whether this increased baroclinity is 
real, or l measure of noise in the data set. (In this case it is the latter.) 

The next set of comparisons involve visual inspections. Analyses of 
the 300- mb geostrophic wind speed are provided in Figure 19. In all three 
satellite data sets, the major jet maximum in the southeast U.S. is placed . 

too far south and west, with an extension back into New Mexico. Also, geo- 
strophic wind calculations in the southwest grid comer are unreliable due 
to a small Coriolis parameter and, more importantly, problems with interpo- 
lation into the wide data gap. 

Differences between the three satellite sounding sets are accentuated 
in these isotach analyses. The SCR 4x7 set failed to define the jet maximum 
propagating down the west side of the east coast trough, although it did show 
stronger wind 3 aver the St. Lawrence River Valley. Overall, the SCR DST iso- 
tach analysis was probably best. It indicated another jet maximum over the 
Carolinas, and produced a better isotach configuration on the west side of 
the trough. In fact, the 45 m s -1 geostrophic wind maximum west of Hudson 
Bay in this analysis may be a feature not resolved in Canadian radiosonde 
height measurements. A jet maximum was present in NMC's 0000 GMT 300 mb wind 
analysis over Trout Lake, Ontario (54°N, 90°W) , which was due mainly to a 
68 n s' 1 wind reported at that station. Considering the general underesti- 
mation of gradients in the satellite-derived analyses, the actual geostrophic 
speeds of this feature may have been even stronger. 

Geostrophic temperature advections at 850 mb are fairly well repre- 
sented in the satellite-derived analyses (Figure 20). In the Eastern region, 
the 6trong cold advection in northern Georgia and adjoining states is captured 
in all three satellite sounding sets. The SCR 4x7 set provides a better indi- 
cation of the maximum south of Lake Ontario than the DST sets, although the 
central value is still weaker than that from the LFM analysis. 












Figure 20. Geostrophic temperature advection 
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Tlic 200 mb goostrophic temperature advection (Figure 21) further 11- 
lustrat.s the gancr.l proble. aatalute .ouadlag. have !„ resolving th . ,„ po . 
pause. The ee.perature gradient at this level has reversed in direction (to. 
that in the nld-ttoposphete. While »»„. general agreenent be found in 
the southeast, the strength of the warm advection in the northe.at (10‘C da,' 1 ) 
and the large are. of cold advection in the central United State, hav. boon 
grossly underestimated in the satellite sounding data. 

For the final analysis comparison, isentropic cross sections were con- • 
structed across the frontal zone situated along the east coast, as illustrated 
n igure 22A. Analysis values were interpolated from the LFM grid points to 
the radiosonde locations shown. Panel B is from the SCR DST analysis and V 
panels C and D are from the bracketing 1200 GMT and 0000 GMT LFM analyses 
respectively. (The satellite analysis grid did not extend far enough east to 
include the Bermuda station, so values at the intersection of the path of the 
cross section with the last row of the analysis grid were plotted for the SCR 
DST analysis set *> Also Shown are analyses of the thermal wind component " " 
norma! to the cross section, built up from a value of zero at 850 mb. Poai- 
tive values point into the page. While the wind maxima in the SCR DST cross 
section are well positioned horizontally, they are about 50 mb higher than - 
those in the LFM cross section. Also, the eastern maximum is about 35X too 
weak in the SCR DST set. The cold dome is much weaker in the satellite analy- 
sis, and its isentropos show almost no slope near the tropopause, a good vis- 
ualization of the deficiencies of the satellite-derived tropopause data. With ' 
^hjogJU^ s ults near the almost cvcrv awtll „. Mw , 

mioter^ there may b ^ome_au opti on in using satel lite sounding heirht data 
above about 300 mb . ’ 

c. LFM initialization with satellite-derived analysis fields 

In the process of defining model initial fields from observations 
several operations affecting the temperature gradient information are per- 
formed. Consider, for example, the steps taken in preparing the LFM initial 
our fields. Rawinsondo measurements of temperature, pressure, moisture and 
wind are made at significant levels. However, only values at the ten manda- 
tory pressure levels from the surface to 100 mb and tropopause data, enter 
the LFM analysis. . Finally, the resulting analysis fields are used as input . 
for the initialization procedure. 
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Figure 22. Analysis cross section comparisons. Isentropes (K), solid and 
thermal wind component <m normal to the cross section «J«ive to 

850 mb. dashed (positive into the page). Panel A - cross action s at 
i p-_-i o _ crp DST analysis cross section at roughly 1800 jMT. 

Panels 2 and D - bracketing 1200 GMT and 0000 GMT LFM.analysis cross sections. 
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The various steps in this process are illustrated with cross sections 
from 1200 GMT 22 February 1976 in Figure 23. (This is the same cross section 
employed in Figure 22.) Only subtle differences can be found between the sig- 
nificant level and mandatory level analyses. A 13X loss of gradient informa- 
tion is indicated by a decrease in the maximum thermal winds in the step from 
mandatory level data to LFM analyses. Model initialization has only a slight 
effect on the analyzed temperature gradients. The dramatic losses shown by 
Horn et al. (1976) also probably occurred in the analysis step, not in the 
initialization step. 

The differences between the satellite-derived analysis fields and their, 
corresponding initial hour fields were quite small, with only a 2X loss in 
height gradient information. The available potential energy decreased AX for 
the LFM fields and only 2Z for the satellite fields due to initialization. \ 
In both the satellite and conventional data cases, mass field changes due to 
the initialization process were small. The smooth vertical nature of the 
satellite-derived data provides little benefit to the initialization proce- 
dure. * 

6. Summary and conclusions 

Detailed comparisons between analyses constructed from satellite sound- 
ings and surface ’data only with conventional analyses from the LFM model, have 
been presented for one period from DST-6. The grid tnesh employed in this 
study has a higher resolution (190.5 km) than those used in most of the pre- 
vious DST-6 impact studies (A00 km) . As stated earlier in the description of 
the experiment, these satellite sounding evaluations addressed questions re- 
lated to the following topics: defining the strengths and weaknesses in sat- 

ellite soundings, investigating the effect of increasing the horizontal reso- 
lution of satellite soundings, assessing the impact of careful data screening 
on the final sa tel lite-derived analyses, and evaluating the effect of model 
initialization on the analyses. The results of this study in regard to these 
topics are summarized in the following discussion. 

Many of this study's findings parallel those from the coarser resolu- 
tion DST-6 impact tests of Dcsnarais et al. (1978) and Halem et al. (19 .’>),• 
and tests with data from earlier satellites, such as Bonner et al. (1976). 

The satellite soundings are able to define the major trough and ridge posi- 
tions quite well, but are ,, conservativc ,, t with temperatures too warm in 
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Fisture 23. A cross-sectional representation of thernal gradient changes during the 
LFM analysis - initialization process at 1200 GMT 22 February 1976. The parameter, 
end station locations are the sane as those in Fig. 21. The upper panels used c g- 
nif leant and mandatory pressure level radiosonde reports. The lower panels are for 
the LFM analysis and initialization fields. 
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troughs and too cold in ridges. In fact, the results presented here indicate 
that the inability of Nimbus-6 satellite soundings to define cold domes out- 
weighs problems in the ridge regimes during DST-6. 

The 4x7 data set presented here has a higher horizontal resolution 
(about a 150 km spacing) than the sounding sets provided for the DST impact 
studies (about 300 km). The higher noise level and resulting higher removal 
rate in the screening procedure of 292 for the ALL 4x7 set compared to 52 for 
the ALL DST set, suggest that both good and poor information from 4x7 data 
blocks were averaged together to yield mediocre yet consistent DST soundings. 
While the manual screening of the ALL 4x7 set produced pronounced improve- 
ment^ the SCR 4x7 set still had a slightly higher noise level and a less 
uniform observation distribution than the SCR DST sounding set. These fac- 
tors are probably responsible for the disappointing SCR 4x7 performance. How- 
ever, the removal of poor 4x7 block data before the averaging to produce DST 
soundings may have yielded better results. Only minor improvements could be 
noted from screening the ALL DST data set. LFM initialization with these 
satellite data sets is feasible, and has only minor effects on thermal gra- 
dients in the step from analysis to initialization for both the satellite and 
conventional data sets. 

The deterioration of the longwave infrared (15 pm) channels of the HIRS 
instrument before the start of DST— 6 had a detrimental effect on soundings 
used in this study and in the DST-6 impact studies. The seven channels lost 
comprise almost half of the temperature sounding' channels available from 
Nimbus-6 and would provide additional tropopause and tropospheric information. 
Petersen and Horn (1977) showed promising results in tracking a closed 500 mb 
low and its associated wind maximum across eastern Canada from an August 
DST-5 sample of Nimbus-6 soundings before the HIRS instrument malfunction. 

The inability of the Nimbus-6 soundings in this DST-6 study to correctly 
position the jet streak in the southeastern United States could be attributed 
to a number of factors, including the unfortunate position of the data gap in 
the region of interest, and the loss of the 15 pm channels. 

More reliable soundings derived from a complete set of infrared and 
microwave radiance measurements have since become available. TIR05-N and 
N0AA-6 satellite soundings were incorporated into the FCGE year data sets* 
and the quality of the soundings has improved to the degree that they are now 
included in the NMC operational data base over oceanic regions. However, 
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results from Phillips et al. (1979), Schlatter (1961) and Schmidt et al. (1981) 
Indicate that satellite versus conventional differences for these never satellites 
exhibit much the same type of structure as the Hinbus-6 differences presented 
here, especially with large differences near the surface and tropopnuse. A re- 
cent study by Schmidt et al. (1981) demonstrates a conservative nature (too want 
in troughs, too cold in ridges). The similarity in the structure of these dif- 
ferences suggests limitations in the basic approach of inferring temperatures 
from radiance measurements, in radiance measurement accuracy, and/or in the 
procedures used to estimate temperatures from radiances. 

While the overall results from this DST-6 case study are somewhat dis- 
appointing, the soundings verr able to define the major trough and ridge posir 
Lions. Their ability to resolve smaller scale features was inconsistent. 

Tiieie was seme indication that satellite soundings could define certain short 
wave features over Canada that were not evident in the LFM analyses. On the 
other hand other short wave features found in the satellite data were incon- 
sistent with conventional data. • . * 

Satellite soundings are currently unable to define temperatbre features, 
particularly important inversions such as the tropopause, with a detail com- 
mensurate with the radiosonde. One should remember however, that a major 
purpose of satellite soundings is to supplement conventional data in data- 
sparse regions, not to replace the current rawinsonde networks. Considerable 
effort should be focused on further defining what satellite soundings can con- 
tribute, and how they can be mixed with data from other sources. 
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Evaluation of TIROS-M and NOAA-6 Satellite Daft 
Comparisons of Colocated Soundings and Analyses for a January C ace 

Brian D. Schmidt, Thomas L. Koehler and Lyle n. Horn 
Department of Meteorology 
University of Wlaconsln-Madlson, HI 53706 


Abstract 

This study uses an early January 1980 synoptic situation In evaluating the 
performance of TIRCS-N and NCAA-6 operational temperature soundings. Visual and 
statistical comparisons of temperature and thickness fields were employed to deter- 
mine the effects of manual screening of the satellite soundings and measure the 
accuracy of both the satellite soundings and analyses derived from them. A com- 
panion study by Dcrber ct al. (1981) used these satellite data as numerical weather 
predict ion' model input. In both studies, comparisons between the performance of 
T1R0S-N and N0AA-6 were emphasised. 

Although the manual screening did not Improve the statistics, removal of the 
poorest soundings produced more consistent analyses. Both satellites were able to 
correctly position the major troughs and ridges. Crsdlcnts were underestimated 
though, with troughs markedly too warm and ridges ollghtly too cold. The poorest 
data occurred near the surface and tropopausc ns reflected by larger standard 
deviations In those layers. As Indicated In both types of comparisons, the perfor- 
mance of TIROS-N was slightly superior to that of NOAA-6. The results of this 
study are encouraging, although more attention should be directed toward correcting 
the problems satellites exhibit near the surface and the tropopausc. 
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l • Introduction 

k major development in the meteorological community during the 1970’m was 
the ndvont of satelllte-dcrtvcd vertical temperature soundings* These soundings 
promised to Improve numerical weather prediction models by providing an adequate 
data base over the large data-aparae regions of the earth* However* the ability of 
satellite soundings to actually Improve prediction models has been Inconclusive* 
Studies conducted at the National Meteorological Center (NMC) Indicated that the 
addition of satellite soundtnga did not Improve numerical model forecasts (Tracton 
and Mcrheraon, 1977 and Dcntaarnla et al., 1978). Similar studies carried out by 
Halea et al. (1978) and Kelly et al. (1978) reported slightly Improved model fore- 
casts when satellite data were Included. The satellite data In these studies were 
mixed with data from conventional sources complicating the evaluation of the *• 
satellite soundings, 

K second approach in evaluating satellite soundings involves comparing 
them against Independent verification data usually based on radiosonde obser- 
vations. Such evaluations are far simpler and less expensive than the numerical 
model Impact studies noted above. For example, Phillips et al. (1979) examined 
the nature of TIROS-N soundings under varying degrees of cloudiness through 
comparison with colocated radiosondes. Another examples of this second approach 
examines the ability of satellite data to define. synoptic features. Horn ot al, 
(1976) found that Nlmbus-5 soundings successfully located the thermal gradient 
beneath an upper tropospheric Jet streak. TIROS-N data was used recently by 
Strelt and Horn (1981) to track the polar and subtropical Jets over the eastern 
Pacific preceding the Wichita Falla tornado outbreak. Other satellite sounding 
studies Involve constructing height analyses. For example, Petersen and Horn 
(1977) tracked a closed 500eb low using Nlmbun-6 data. More recently, Koehler 
(1981) was able to delineate atmospheric trough and ridge position* over North 
America from Nlrabus-6 height analyses constructed on a fine mesh model grid. 

This study continues the approach of not mixing satellite and radlononda 
data, noth colocation comparisons similar to Thill <s ct al. (1979), and com- 
parisons of synoptic features slmtlar to Koehler (1 • ’l) are performed* Unlike 
previous studies, two satellite data sets (TIROS-N and NOAA-6) will be examined. 

A case study of an early January 19«o synoptic situation occurring in the United 
Staten (t)-S.) will be used. The choice of this area allows the use of the 
familiar United-area Fine Mesh (l.FM) model analysis fields for verification. 
Uerlved entirely from conventional observations, the LFM analysis field* provide 
nn Independent standard of comparison. 
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With TIROS-N and NGAA-6 In operation at that time, a unique opportunity 
was provided to evaluate and compare the performances of the two satellites. The 
evaluation Is pursued In two distinct ways. This study examines the quality of 
both the actual satellite temperature soundings and analysis fields based on 
them. In a second study by Dcrber ct al. (1981), the resulting T1R0S-H and 
NOAA-6 analyses are used to initialize and run a numerical weather prediction 
model to further evaluate the satellite soundings. Thus these two studies employ 
some aspects of the two basic approaches for satellite sounding evaluation. Inde- 
pendent satellite sounding comparisons against conventional data and assessment 
of satellite data as numerical model input* In both studies comparisons between 
TIROS-N and NOAA-6 will be emphasized. 


2 • The region and period of study 

A subset of NMC*s LFM model output grid was chosen as the region of study 
as shown In Figure l. This polar stercographlc grid located over North America 
has a grid spacing of 190.5 kilometers true at 60 # N and is fine enough to delin- 
eate small scale synoptic features revealed by the satellite soundings. 

Reasonably reliable analyses derived entirely from conventional rewinsonde data 
are available on this grid from NMC. In this study they served two purposes! 

1) they provided the verification data set against which satellite sounding data, 
were compared and 2) they were used ns a first guess In an Interpolation scheme 
employed In the horizontal analyses of the satellite data. LFM analyoes were 
extracted from the entire domain shown In Figure 1. However, the satellite- 
derived analyses were constructed only on the sub-region outlined. This analysis 
region was chosen to facilitate Interpolation to the model grid used in the 
companion study by Derber et al. (1981). 

In order to use and evaluate the satellite data, the period of study had 
to satisfy several requirements. The ability of satellite soundings to define 
strong temperature gradients could lest be evaluated In a synoptic situation that 
exhibited Intense baroclinity. A wintertime pattern satisfies this requirement. 
Furthermore, two synoptlcnlly active regions, e.g. well developed troughs, within 
the analysis region would allow for more meaningful evaluation of the horizontal* 
error structure. Also, a sustained period of nearly complete data sets froe V 
TIROS-N and NOAA-6 was necessary for intcr-satelllte comparisons and inter- 
polation to synoptic times. 
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Figure 1« The LFM grid and the analysis subset grid (dark outline) for this 
study. 

i 

A period between 0000 GMT 5 January and 0000 GMT 8 January was found to 
satisfy these requirements. Figure 2 depicts the 500 mb and 1000 nb height 
fields for the period 0000 GMT 6 January through 0000 GMT 7 January 1980. At 
0000 GMT 6 January, a long wave trough was located over the eastern portion of 
the U.S., while a short wave trough was located over the Pacific Northwest* An 
Intense Jet streak was associated with this vigorous short wave as it rapidly 
propagated eastward. Reasonably complete TIR0S-N and N0AA-6 satellite data oats 
were available from. the Environmental Data and Information Service (EDIS) t nndLFH 






ORIGINAL PAGE IS 
°F POOR QUALITY 
















-56- 


save tapes were furnished by NMC for this period* These data sets allowed for 
Interpolation of the asynoptlc satellite data to the synoptic tines, 0000 CUT and 
1200 COT 6 January and 0000 GOT 7 January. 


3. The satellite pounding data 

TIR0S-N and N0AA-6 operational satellite soundings were used In this study. 
The character of the study was Influenced by the satellite orbits. Instrument 
design and retrieval procedures. Features pertinent to this study are provided In 
the following discussion. A more detailed description of the orbital and Instru- 
mental characteristics can be found in Kldwell (1979). The operatl >nal satellite 
soundings produced by the National Environmental Satellite Service (NESS) 
contained three major types of soundings, differing In both the raw data and the 
processing procedures. These soundings produced reasonably complete data coverage 
across the analysis region# 

The TIR0S-N and NOAA-6 satellites were launched on 13 October 1978 and 27 
June 1979 , respectively# They were placed in nearly sun-synchronous orbits with 
southbound equatorial crossings at approximately 0300 (TIROS-N) and 0730 (NQAA-6) 
Local Solar Tice (LST) • With northbound equatorial crossings occurring twelve 
hours later f the analysis region was covered twice each day by each satellite. 

For example* Figure 3 displays sounding locations obtained from TIROS-N on 6 . 
January 1980- The orbital period of approximately 102 minutes resulted in an 
observation time variation of nearly ten hours from east to west across the analy- 
sis region- This asynoptlc nature of the data made a time interpolation necessary 
to produce analyses at standard times to facilitate verification# 

Both satellites carried similar versions of the TTIR0S Operational Vertical 
Sounder (TOVS). The TOVS consists of three instruments* the High resolution 
Infrared Radiation Sounder (HIRS-2) , the Microwave Sounding Unit (HSU) and the 
Stratospheric founding Unit (SSU)- Vertical temperature profiles produced by 
NESS were derived from the radiance measurements provided by these instruments# 
Clear, partly cloudy and cloudy soundings were produced using procedures described 
in Smith et al. (1979). Clear soundingn were derived from a combination of all 
available HIRS and KSU channels. A special method described by Smith and Woolf 
(1976) which compensates to some degree for cloud contamination of the Infrared 
(HIRS) channels was used for the partly cloudy soundings. Cloudy soundings were 
derived from all stratospheric channels* but only the microwave channels in the 
troposphere. 
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Pig. 3. TIROS-N sounding locations from descending orbits at approximately 
0300 Local Solar Time (LST) on 6 January 1980. Soundings are divided into three 
types: clear (0), partly cloudy (a) and cloudy (+). 

At the time of this case study (January I960), the acceptance criteria for- 
the partly cloudy retrieval technique was such that too few partly cloudy soun^ 
dings were produced (Smith. 1981). Also, the automated quality control procedure 
used during the data processing was too liberal, making manual screening of the 
satellite soundings important. Previous work by Uorn et al. (1976), Blechman and 
Horn (1981) and Paulson and Horn (1981) has shown that such manual screening can 
produce data sets superior to unscreened satellite data sets. 


4* Data preparation and screening 

The original format of the operational satellite soundings received from 
EOIS was in terms of layer mean temperatures, along with temperatures and 
pressures at the surface and tropopause. This study's screening and analysis pro- 
cedures required that the data be in terms of lsobaric level temperatures and 
thicknesses relative to lOOOmb. Using a method described by Polger (1978), the 
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level temperatures were determined by a linear Interpolation between the. layer 
mean temperatures, surface temperature and tropopause data* Level temperatures 
below the reported surface were estimated assuming a lapse rate of one-half the 
dry adiabatic lapse rate* Thicknesses relative to lOOOab were hydrostatically cal* 
culated using the original layer mean temperatures where possible, and the layer 
temperatures estimated from the level values below the surface and near the trope- 
pause* 

For screening and analysis purposes, the data were divided Into five sets 
per satellite, each set consisting of soundings from either ascending or descend- 
ing orbits at the same subsatellite LST* A sample of one of these data oeto is 
given in Figure 3. The arduous process of manually screening the ten data seta 
was undertaken by intercomparing horizontal hand analyses, vertical cross auc- 
tions, cloudiness and type of sounding* Initially, questionable soundings were 
flagged in hand analyses of temperatures at seven mandatory levels op to ICOmb and 
also for cumulative thicknesses (i#e* 1000-850mb, 1000-700nb, etc*)* Vertical 
cross sections of potential temperatures and thermal winds relative to 830®b were 
constructed through the flagged data points using a technique developed by v 
Whittaker and Petersen (1977). When the horizontal analyses and vertical cross 
sections failed to clearly identify inconsistent soundings, other factors such as 
cloudiness and type of sounding were considered. When this occurs, the cloudy 
(microwave) soundings and soundings with higher percentages of cloudiness were 
deemed inconsistent* To reduce the effect of personal bias, agreement between the 
two people involved in the screening procedure was necessary* During this pro- 
cess, direct comparisons between satellite 'and conventional analyses wars avoided 
to insure independent judgements* It was found In the screening procedure that 
Inconsistent soundings which were too cold (warm) in the lower layers were too vara 
(cold) in the upper layers. In regions of overlapping passes, sounding pairs lata 
than sixty kilometers apart could produce instabilities in the analysis routine. 
Special attention was given these pairs with the less consistent sounding being 
removed from all data sets. 

From these screening procedures, ten data sets for each satellite emerged, 
five screened (i.e. obviously inconsistent soundings removed) and five unscreened 
(i.e. obviously inconsistent soundings retained). For TIROS-H, the unscreened data 
sets contained a total of 2095 soundings. Of these, 6.4 percent were removed to 
create the screened data sets consisting of 1962 soundings. Likewise, the ftGAA— 6 
unscreened data sets contained 2257 soundings with 6.4 percent removed to form 
screened data sets consisting of 2113 soundings* 
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5* The asynoptlc problem 

A major difficulty encountered In studied coopering satellite soundings 
with conventions! data is the time variation between the two types of data* Con- 
ventional observations and the analyses based on then are synoptic in nature* 
However, satellite soundings are made in continuous orbital swaths at the sane LST, 
resulting in tine differences of 102 minutes between consecutive passes* To cover 
the Horth American region, at least six passes are roqulred, which leads to time 
differences of nearly ten hours between the easternmost and westernmost soundings* 
Since this type of study requires comparisons of the asynoptlc satellite data with 
an independent verification set, both fields oust be valid at the same time to pro- 
duce truly meaningful comparisons* To place both data sets in the same tine frame 
for comparison, one or the other must be modified* If the comparisons are per- 
formed at the satellite observations, a scries of bracketing conventional standards 
of comparison should be Interpolated to the satellite times* This allows statisti- 
cal comparisons which identify actual sounding error characteristics* On the other 
hand, if the comparisons are performed at the synoptic times of conventional analy- 
ses, a series of bracketing satellite f lelds shouli be Interpolated to synoptic 
time* This provides an opportunity to evaluate the performance of satellite data 
in defining the structure of weather patterns. However, time interpolations nay 
provide an inadequate representation of the actual temporal changes in the 
atmosphere* 

To facilitate an understanding of this complex topic, a detailed description 
of the time-related aspects of this study will be presented* An appreciation of 
this topic is essential to better understand the results given in later sections* 

Conventional NMC analyses on the LFM grid based on synoptic radiosonde 
observations provide the standard of comparison for the satellite sounding eva- 
luations* Assessment of satellite data is performed using both types of com- 
parisons noted earlier: 1) statistical comparisons done at satellite sounding 

locations by interpolating the gridded LFM data in time and space to the sounding 
locations (colocated comparisons) and 2) comparisons of synoptic features based oa 
satellite and LFM analyses at synoptic times. 

Figure 4 illustrates the time dependency of the data sets employed. Both 
the discrete synoptic times of the LFM analyses and the continuous periods of the 
satellite ooundihgs are presented. For the colocated comparisons, the seven LFM 
analysis periods shown schematically in Figure 4 produce verification values for 
the four shaded TIROS-N and NOAA-6 data sets. Both morning and afternoon coloca- 
tions were performed to prevent a diurnal bias. 
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Figure 4. The satellite data sets and LFM analyses used in this study. 
Comparisons at sounding locations were made for the shaded satellite sets. 
Satellite-derived analyses were interpolated to the times shown at the 
bottom for synoptic feature comparisons. 



Figure 5. Satellite sounding- tine analysis for the 6 January 1980 
N0AA-6 corning soundings, with observations before 1200 GMT shaded. 
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Verlf.lcatlon values for each, level at every satellite sounding location were 
determined aa follows. First, the satellite sounding time was used to identify 
four bracketing LFM times. A spatial interpolation for each of these four analyses 
was then performed to the satellite sounding location using overlapping quadratic 
polynomials in two dimensions. This was followed by an overlapping quadratic poly- 
nomial interpolation in tine to the sounding time, producing the colocated value. 
This process is clarified by considering a complete satellite period as shown io 
Figure 5, which presents a time analysis of the 6 January 1980 NCAA-6 morning soun- 
dings. LFM periods B, C, D and B (Figure 4) were used in Interpolations for all 
observations in the shaded area (i.e. before 1200 GMT), while LFM periods C, D, E 
and F were employed in interpolations in the unshaded area (i.e. after 1200 CMT).. 
Thus, the production of colocated verification values for all four TIE0S-M and 
NQAA-6 periods required seven LFM analysis sets. 

Conversely, for comparisons at synoptic times, five satellite periods werq 
employed to construct the three satellite synoptic analysts shown in Figure 4. 
First, satellite analyses at LFM grid points were constructed for each of the flva 
TIR0S-H and NQAA-6 data sets. Ignoring time variations in the observations. 

Separate time analyses were also constructed, such as that shown in Figure 5. To 
deternlno the ’synoptic’ satellite value at an LFM grid point, the synoptic tio* 
was compared to the values of the five time analyses, the bracketing times for the 
polynonial interpolation identified, and the value calculated. This process was 
repeated at each level for every grid point, yielding the fiosl satellite synoptic 
fields for the three times shown in Figure 4. In any event, oore data sets were 
'involved in generating compatible comparisons than were actually evaluated in the 
result sections. 


6. Thickness analysis procedure 

Satellite thickness analyses were produced on the LFM grid using an optimum 
Interpolation method with normalized weights. T‘ is method, developed by Koehler 
(1979) based on Candin (19*3), w«s specifically designed to handle the abrupt 
changes in coverage, such as gaps between satellite passes. Since the inter- 
polation weights are normalized, the analysis can be performed either with or 
without first guess fields. Persistence first guess fields provided by the LFM 
thickness analyses approximately twelve hours preceding the average sounding ties 
were used in the satellite thickness analyses. Thickness analyses were produced at 
nine vertical levels with identical weights used at each level. The satellite time 
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analyse** required for the Interpolation to synoptic times* were const leered 
without a first guess using the same weights employed in the thickness analyses* 

The combined filter described by Cerrlty and Newell (1976) was applied to the 
final* tine interpolated, synoptic satellite thickness analyses* 

7* Results 

The discussion of the results is divided into two sections: 1) colocated 

comparisons and 2) the synoptic pattern analyses. Because this study examines the 
performance of two satellites over a two day period* a thorough discussion of the 
results could become tedious. Thus* while detailed results are presented* only 
those features deemed most significant will be discussed. 

a) Colocated comparisons 

As previously mentioned* this approach compares satellite layer mean tem- 
peratures with conventional gridded data interpolated to the satellite observation 
locations. This comparison provides an estimate of the accuracy of the soundings* 

Table t shows TIROS-N and NOAA-6 sounding statistic* for the screened and 
unscreened data sets. Since there are only very minor differences between tho two 
sets of data* it is apparent that manual screening did not improve the overall sta- 
tistical results. However* it must be emphasized thai the screening procedure does 
remove the poorest soundings* thus providing more consistent satellite sounding 
sets for analysis purposes. For this reason* all future comparisons will use only 
the screened satellite data sets. 

Figure 6 shows a comparison of TIROS-N and NOAA-6 sound .ngs over both the 

whole region cf study and the U.S. subregion. The data in this figure ara a 

composite of all four tine periods and all three soundings types. Excluding the 

lower two layers* the graphs exhibit strikingly similar patterns for both fatal- 

lites. The larger bias and standard deviation values in the 200-300nb region 
reflects the difficulty in defining the tropopauae using satellite soundings. 

The U.S. region exhibits a much. more pronounced pattern* with larger negative 
biases in the upper troposphere. Several factors may contribute to this pattern* 
The U.S. region with its dense radiosonde network* provides a more stringent 
verification standard than that available in the entire region. Additionally* 
the U.S. region la primarily in mid-latitudes where January baroclinlty is 
strong, while the whole region lrcludes areas north and south of this strongly 
baroclinic region. In these weaker barocllnic areas there is. less likelihood of 
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Figure 6. Colocated layer mean temperature comparisons (*C) between 
TIROS-K/NOAA-6 values and LFM analyst* values (interpolated in space 
and tine to the satellite sounding locations). 
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Table Is Layer mean temperature comparisons (*C) between screened (SCR) and 
unscreened (ALL) satellite soundings and analysis values Interpolated to obser- 
vation locations* 


LATER 



1000- 

850 

850- 

700 

700- 

500 

500- 

400 

BIAS 

ALL NOAA-6 

1.08 

.83 

.04 

-.36 

SCR NOAA-6 

1.17 

.89 

.08 

-.34 

ALL TIROS-N 

-.12 

.24 

.08 

.06 

SCR TIROS-N 

-.01 

.30 

.12 

.07 

ST. DEV. 

ALL NOAA-6 

4.51 

2.68 

2.29 

2.20 

SCR NOAA-6 

4.47 

2.64 

2.28 

2.18 

ALL TIROS-N 

4.16 

2.40 

2.04 

2.08 

SCR TIROS-N 

4.12 

2.33 

2.02 

2.09 
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400- 

300 

300- 

250 

250- 

200 

200- 

150 

150- 

100 

-.10 

1.59 

1.35 

.13 

-.14 

-.09 

1.58 

1.29 

.09 

-.15 

.32 

1.74 

1.42 

.21 

.02 

.32 

1.74 

1.41 

.20 

.03 

2.46 

2.69 

3.45 

2.89 

2.17 

2.48 

2.70 

3.45 

2.88 

2.17 

2.45 

2.51 

2.81 

2.22 

1.94 

2.46 

2.52 

2.79 

2.20 

1.92 


large differences between the satellite soundings and radiosonde based verifica- 
tion data. Thus with the larger sample, the differences should be reduced* 
Finally, within the larger north-south extent of the whole region, there Is a 
greater variation in the height of the tropopause. Thus the distribution of dif- 
ferences is spread over several layers, which reduces the magnitude of the dif- 
ferences in the upper troposphere* . . 

Figure 7 presents the TIROS-N and NOAA-6 soundings subdivided into the 
three sounding types: clear (CLR) , partly cloudy (P CLDY) end cloudy (CLDT). 

This type of figure, similar to one used in Phillips et al* (1979), reveal* the 
characteristics of the various sounding types. 

In the Phillips study, TIROS-N soundings from various regions of the globe 
during April 1979 were compared with colocated radiosondes* Bias and RMS values 
were calculated* In this study, the standard deviation is used, rather than. the 
RMS. The standard deviation provides a better measure of the "noise** since the 
bias has been removed. 

A comparison between the bias values obtained for TIROS-N In this study 
(based on the entire region) with those obtained by Phillips et al. (1979) for 
the North American region shows similar values with largest biases amounting* to 
about 2*C. While the largest biases in this January study occur at about 200- 
300mb, the largest values in the Phillips et al. study occur between 150-250ab. 
This probably reflects the higher altitude of the tropopause in April. Overall, 
the bias values obtained by Phillips et al. for North America are not greatly 
different from those in this study based on the entire region. However, biases 









based on the U.S. subregion show considerably larger values than obtained by 
Phillips et al.. The largest standard deviation near the tropopause obtained In 
this study (based on the entire region) of 2.95*C yields an RMS of 3.3*C when 
combined with the bias. This compares well with an RMS value near the tropopause 
of about 3.4*C in the Phillips et al. study. On the other hand, an RMS of 4.8*C 
la reached for the U.S. subregion of this study. However, the standard devia- 
tions for both TIROS-N and NOAA-6 in the U.S. subregion tend to be somewhat 
smaller than those In the whole region. This stands In contrast to the larger 
bias values In the U.S# subregion. 

As noted in the Introduction, the major aspect of this study Is the com- 
parison of the relative performance of two satellites, TIROS-N and NOAA-6. 

Figure 8 provides a direct comparison of TIROS-N and NOAA— 6 for the clear and ' 
cloudy soundings over both the entire region and the U.S. subregion, ihe partly 
clcudy soundings are not presented since they comprise such a small percentage of 
the sample. As shown in Figure 8, the TIROS-N and NOAA-6 bias and ctandard de- 
viation curves are quite similar above the lowest layers* Thi9 similarity is not 
surprising since pearly identical instruments are aboard both satellites* .How- 
ever, the standard deviations show NOAA-6 to be slightly noisier throughout the 
soundings a 

The most significant difference in the curves occurs in the bias for the 
1000-850mb layer where TIROS-N soundings are noticeably colder. In an attempt to 
isolate the reason for these low level differences, histograms were prepared 
showing the distribution of the TIROS-N and NOAA-6 colocated temperature dif- 
ferences over the eastern' U.S* (Figure 9). This region of the U«S* was used to 
lessen the problems of topography. The NOAA— 6 clear cases (Figure 9a) shows a 
distinct shift toward warmer values (l.e. NOAA-6 warmer than the radiosonde)* In 
contrast, the TIROS-N clear cases (Figure 9b) are okewed toward values colder 
than the radiosonde. The cloudy comparisons provided in Figures 5c and 9d show a 
strong skewness toward warmer temperatures at low levels for NOAA-6 while the 
TIROS-N data are more normally distributed. Attempts were made to explain dif- 
fering distributions shown in Figure 9 without success. Perhaps future studies 
should further Investigate this tendency. 

b) Height and thickness synoptic analysis comparisons 

Figures 10 and 11 depict the TIROS-N, NOAA-6 and LFM 500nb height 
and geoatrophic isotach analyses for all three time periods. Both satellites 
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Figure 11. Same as Fig. 10 except for NOAA and LFM. 
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define the trough and ridge positions relatively well, althou B h none of the gra- 
d lints ore weaker than In the LFH analyaea. Thla la particularly evident at 0000 
CHT 7 January, when both aatellltea fall to capture the 60 m a Jet naxlnua 
associated with the trough located In the Northern Plain, states. Weaker 
satellite gradients have been noted In early studies (e.g. see Koehler, 1981 or 
Oesaarals at al.. 1970). Also, at 1200 CHT 6 January both satellites show a ten- 
dency toward a double wave In the western trough. This tendency Is not present 
in the LFH analysis. A further discrepancy occurs over Mexico, where apparently 
anomalous troughs are found in several of the satellite analyses. These trough. 

.re located In data gaps between satellite passe, and can probably be attributed 
to the extrapolation of gradients Into the data gap. by the analysis method. 

In summary, an examination of the 500mb height, and geo.trophlc wind, 
shown In Figures 10 and 11 are relatively encouraging despite the discrepancies . 
noted above. The quality of these satellite analyse, will be further revealed In 
numerical forecast experiments conducted by Berber et al. (1981). 

Since the colocation results showed relatively poor value, at low level, 
and near the tropopause, an attempt was made to decrease the effect, of these • 
layers In further evaluations. This was done by calculating thicknesses over •• 
three tropospheric layers (i000-700mb, 700-300mb, 300-l00mb). The discussion 

that follows pertains to these three layers. 

The satellite minus LFM 1000-700mb thickness differences shown in Figure L 
illustrate satellite sounding problems near the surface. Both positive and nega- 
tive centers arc present. Over southern Canada and the O.S., both satellites show 
positive centers (l.e. satellite thicknesses too high) propagating eastward and 
southward with the major trough. This is consistent with the 500mb trough 
(Figure. 10 and 11) being coo weak. On the other hand, the negative anomaly (l.e. 
satellite thicknesses too low) remains more stationary, primarily over the moun- 
tainous western region of the U.S. This may Indicate a recognised difficulty In 
the retrieval procedures over high terrain, or It may reflect inadequacies In the 

extrapolation to lOOOnb In this study. 

Two other large anomalies exist in northern Canada. One Is a large posi- 
tive anomaly located over the North Atlantic coastal region near Labrador. Since 
It is located near the center of the North American branch of the polar vortex and 
occurs In both satellites it may simply represent another example of trough, being 
too warm In satellite analysis fields. However, since thla are. Is nearly devoid 
of conventional data, It la possible that the anomaly reflecta a deficiency In the 
verification data. The other anomaly exists over, mountainous veatern Caned, and 
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Figure 12. Satellite-derived minus LFM thickness differences (data); 
for the l000-700mb layer from the 0000 GMT 6 to 0000 GMT 7 January 
1980 periods. 
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ia also positive. This feature could be another example of satellite sounding 
problems over high terrain. However, since It weakens with time it may also 
reflect the propagation of the trough from this region. 

The largest departures are in the NOAA-6 data, as compared to the TIROS-H 
data. Since both data sets are Interpolated to the same time, the greater NQAA-6 
anomalies may reflect problems with either its Instrumentation or sounding 
retrieval procedure. Thle la also supported by the slightly better performance of 
TIR0S-N in the colocation comparisons* 

Results for the 700-300ab layer should be more promising, considering the 
smaller colocated differences found in that layer. Figures 13 and 14 chow the 
TIROS-N, NOAA-6 and LFH 700-300mb thickness analyses for the three time periods. 
Also shown are lsotachs (daBhed lines) of the 700-300nb thermal wind. The trough* 
In the thickness field are located west and north of their position In the 500mb 
analyses, which would be expected In a developing barocllnlc wave. Many of the 
variations present In the 500mb fields are removed by Isolating this layer from 
the lower levels. For example, both the TIROS-N and NOAA-6 700-300mb thickness 
patterns give little or no Indication of the tendency toward a "double wave" pat- 
tern In the 1200 GMT 6 January 500mb field. Also, the eastern Pacific closed high 
and low In the 500ob field apparent In all three periods Is weakened In the 700- 
300ob thickness field* 

Both satellites show general agreement with the LFM analysis trough and ‘ 
ridge positions, although the satellite features are shifted 3-5* to the west. 

The Intensity of the troughs and ridges are underestimated with the troughs also 
showing less curvature. As usual, there la a loss of gradient Information In both 
sets of satellite data. 

However, TIROS-N does a better job In defining the 700-300mb gradients than 
NOAA-6. This Is particularly evident In the western trough at 1200 GMT 6 January 
and 0000 GMT 7 January when TIROS-N better defines the 40 m s -1 thermal wind jet 
and related thlcknesB gradient apparent In the LFM analyses. Neither satellite 
does very well In delineating the thermal wind jet associated with the eastern 
trough, although it Is unclear whether this Is totally a satellite problem or 
8 Imply a lack of conventional data over the Atlantic. 

Examination of the statistics In Table 2, which Includes SI scores, also 
reveals NOAA-6' a deficiency In defining gradients. 1 The NOAA-6 SI scores are con- 
sistently larger than those for TIROS-N for both the entire region and the D.S. 


l) 8C0 ™ thfi accuracy of the intensity and positioning of gradients 

centation 0 ° * w th snaller SI scores indicating more accurate gradient repre- 
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Figure 13. T1ROS-N and LFM 700-300ob thickness (solid, dam) and 
thermal wind isotach analyses (dashed, ms*) for 0000 GMT 6 to 
0000 GMT 7 January 1980. 
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Table 2: fIROS-N (T) and NGAA-6 (N) vs* LFM layer mean temperature analysis 

stAtlstlcBi Bias and standard deviation values are in Celsius degrees* while 
the Si scores are dimensionless* 


0000 CKT 6 JANUARY 1200 GMT 6 JANUARY 0000 GMT 7 JANUARY 


ENTIRE U.S. ENTIRE U.S. ENTIRE 0.S* 

REGION SUBREGION REGION SUBREGION REGION SUBREGION 



T 

N 

T 

N 

T 

N 

T 

N 

T 

N 

T 

N 

BUS 













1000-700HB 

-.20 

.83 

-.41 

.81 

-.08 

.29 

-.87 

.45 

-.29 

.74 

-.50 

1.17 

700-300MB 

.09 

.17 

-.77 

-.86 

.15 

-.15 

.92 

-1.38 

•12 

• 10 

i 

. 

<yv 

1 

. 

03 

a* 

300-100MB 

.77 

.47 

1.38 

1.54 

.65 

.35 

1.74 

1.31 

.77 

.37 

1.52 

.87 

ST. DEV. 













1000-700MB 

2.20 

2.82 

1.70 

2.41 

2.59 

2.59 

2.59 

2.50 

2.50 

2.69 

2.23 

2.85 

700-300MB 

1.36 

1.64 

1.17 

1.65 

1.39 

1.45 

. 1.48 

1.54 

1.30 

1.58 

1.10 

1.35 

300-100MB 

1.31 

1.84 

1.06 

1.74 

1.52 

1.91 

1.10 

1.70 

1.48 

1.98 

.91 

1.47 

SI 








* 





1000-700KB 

51.0 

60.2 

36.9 

47.5 

50.6 

51.6 

42.6 

45.0 

49.4 

54.1 

36.9 

45.3 

700-300MB 

42.2 

45.9 

32.6 

39.4 

43.4 

45.1 

35.2 

41.4 

38.9 

42.8 

31.0 

41.0 

300-100MB 

63.5 

74.2 

45.0 

68.7 

69.5 

81.7 

52.2 

79.4 

65.9 

73.8 

47.9 

63.3 
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subreglon. The SI scores for both satellites arc smaller In the U.S* subregion 
then In the whole region* One poosible explanation is that the lack, of radiosonde 
data in the data sparse areas of the entire region prevents an accurate definition 
of the verification field gradients* The standard deviations in Table 2 indicate 
greater values for NOAA-6* This agrees with the colocation results* 

Figure 15 ohows the 700-300ab thickness difference fields (i.c. satellite 
sinus LFM analysis)* As in the 1000-700mb thickness differences* the positive 
anomalies* indicative of too high satellite thicknesses* follow the troughs 
throughout the period in both satellites* This is true for both the western and 
eastern troughs* However* unlike the 1000-700mb thickness difference fields* 
negative anomalies indicative of too low satellite thicknesses* are found in cost 
of the ridges. This is apparent in the ridge between the eastern and western 
troughs throughout the three else periods. Negative anomalies also tend to remain 
stationary over the western O.S* mountains* This possibly reflects the effects at 
higher levels of sounding problems over high terrain. The 70Q-300mb thickness 
difference results support the premise that NOAA— 6 is poorer in defining gra- 
dients. For example* NOAA-6 exhibits a large negative anomaly over Colorado 
at 1200 CKT 6 January and a large positive anomaly over the Northern Plains. at . 
0000 GMT 7 January. 

The 300-100nb analysis results are only given in statistical fora. These 
results are poorer. than those in the 700-300mb layer for both satellites with con- 
sistently larger SI scores and biases. This is not surprising, since the 300— 
lOOsb layer contains the large tropopause errors found in the colocated results. 

8. Conclusions 

Evaluation of TIROS-N and NOAA-6 soundings was performed using both colo- 
cated and analysis comparisons between satellite and conventional data for an 
early January 1980 synoptic situation. Although paralleling similar studies 
(Desmarals et al., 1978; Koehler* 1981; Schlatter, 1981), this study is unique in 
that it simultaneously evaluates the performance of two satellites. 

Extensive manual screening was performed on the original soundings sets. 
Although the colocated statistics showed little or no improvement* removal of the 
poorest soundings provided more consistent analyses. 

The TIROS-N NOAA-6 colocated statistics exhibit similar characteristics 
above the lover layers (l.e. 1000-700ab) , with the best results in the middle tro- 
posphere. The large differences near the tropopause are indicative of poor ver- 
tical resolution in satellites, while those near the surface are probably due to 
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retrlcvisl problems over Itt^h terrain and poor vertical reuclutlon. The TIROS-N 
atatlatlca for the entire region are similar with those obtained by Phillips et al* 
(1979)* However t the TIROS-N and KOAA-6 biases in the U.S. subregion are larger 
than those of the entire region. 

TIROS-N and KOAA-6 analyses ehow general agreement with the LFM SOOab height 
and 700-300ab thickness analyses, especially in locating trough and ridge posi- 
tions. However, some differences were noted. For example, several of the 
variations found in the 500mb satellite fields were not present in the 700-300ab 
thickness fields. Tht large differences found 1 q the layers below 700ab help 
explain the improvement* The 700-300tab thickness differences exhibited positive 
anomalies (coo warm) in troughs and negative anomalies (too cold) in ridges with 
the positive anomalies especially moving with the synoptic situation. These anoma- 
lies Indicate that satellites underestimate gradients and that the underestimation 
is correlated to the synoptic pattern. 

Overall, TIROS-N performs slightly better than NOAA-6 in this study. The 
smaller TIROS-N colocated standard deviations Indicate less noise in the soundings* 
Visual comparisons of the 500ab height and 700-300mb thickness analyses reveal more 
accurate TIROS-N gradients. This is supported by the smaller TIROS-N anomalies in 
both 1000-7 OOab and 700-300mb thickness difference fields. Additionally, smaller 
TIROS-N SI scores quantitatively show wore accurate gradient intensity and posi- 
tioning. 

In summary, the rather accurate description of synoptic features by TIROS-N 
and KOAA-6 Is encouraging, although it is evident that the problems near the sur- 
face and tropopauae have not been corrected. Improvement in these areas should be 
given greater priority* This study* a results show TIROS-N soundings to be of 
slightly better quality than KOAA-6 soundings* Whether this result is true in 
general or is simply an anomaly in this one case study remains to be seen. 
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A Numerical Evaluation of TIROS-N and NOAA-6 
Analyses In a High Resolution Limited Area Model 


John C. Berber, Thomas L. Koehler and Lyle H. Horn 
Department of Meteorology 
University of Wisconsin-Madison, WI 53706 


Abstract 

Vertical temperature profiles derived from TIROS— N and NQAA— 6 radiance 
measurements are used to create separate analyses for the period 0000 GKT 6 
January to 0000 GMT 7 January 1980. The 0000 GMT 6 January satellite analyses 
and a conventional analysis are used to Initialize and run the University of 
Wisconsin’s version of the Australian Region Primitive Equations model. Forecasts 
based on conventional analyses are used to evaluate the forecasts based only on 
satellite upper air data. 

The forecasts based only on TIROS-N or NOAA-6 data did reasonably well in 
locating the main trough and ridge positions. The satellite initial analyses 
and forecasts revealed errors correlated to the synoptic situation. The trough 
in both TIROS-N and NOAA-6 forecasts which was initially too warm remained too 
warm as It propagated eastward during the forecast period. Thus, it is unlikely 
that the operational satellite data will improve forecasts in a data dense 
region. However, in regions of poor data coverage, the satellite data -should 
have a beneficial effect on numerical forecasts. 
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1 . Int r oduct Ion 

Since the advent of numerical weather prediction, the development and uso 
of numerical models have been hampered by the lack of complete and accurate data 
seta, especially in oceanic regic-s. Recent advances in remote sensing have 
established the possibility of using satellite derived temperature profiles to 
Improve the data base. However, in several recent numerical experiments, the 
results have failed to show any consistent positive impact when satellite tem- 
perature profiles were Included in the model initial fields* These results have 
raised questions concerning the usefulness of the satellite data in numerical 
weather prediction. Since many other factors can influence forecast quality, 
this conclusion may be faulty. This study, along with a companion study by 
Schmidt et al. (1981), attempts to determine whether vertical temperature pro- 
files derived from TIROS-N and NOAA-6 radiance measurements contain useful 
meteorological information. In this paper, the evaluation is performed by com- 
paring numerical forecasts made using only sounding data from each satellite 
with a control forecast baaed on only conventional data. 

In previous studies the quality of satellite data has been analyzed using 
two distinct methods: direct comparisons of satellite data with conventional 

data and numerical model forecasts. The direct comparison can either be done 
using colocated soundings or analysis comparisons. Both direct comparison 
methods are used by Schmidt et al. (1981) on the same data set employed in this 
study. The colocated sounding technique compares Individual satellite profiles 
to nearby radiosonde profiles or analyses interpolated to the sounding loca- 
tions. Phillips (1979) and Schlatter (1981) have recently compiled colocated 
statistics for the TIROS-N satellite. In addition, Schlatter (1981) also 
employed analysis comparisons. Ralem et al. (1978) and Tracton et al. (1980) 
combined Nimbus-6 satellite data with conventional data to produce horizontal 
analyses for use in their impact studies. However, the combination of the 
satellite data with conventional data makes satellite data evaluation difficult. 
Peterson and Horn (1977) used only Nimbus-6 data to track a 500mb low. Koehler 
(1981) was able to define major synoptic features using only Nimbus-6 data. 

Using cross-sectional analyses created from Nimbus-5 data, Horn et al. (197.6) 
were able to describe the thermal gradient beneath an intense upper tropospheric 
Jet streak. Blcchman and Horn (1981) used higher resolution Nimbus-6 data to 
better locate a Jet streak in a summertime situation. Recently, using both 
horizontal and. croBs-sectional TIROS-N analyses, Streit and Horn (1981) were 
able to track the’ polar and subtropical jets, which later influenced the Ulchita 
Falls tornado outbreak. 


1 
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Tho inclusion of satellite data In numerical weather prediction asdels 
provides another method of satellite data evaluation. Nearly all of the recent 
studies of this type have concentrated on employing a combination of temperature 
data derived from radiance measurements from Nimbus-6 and conventional upper air 
data. Ghll et al. (1979) found a slight positive Impact with inclusion of 
satellite data, while Tracton et al. (1980) found no consistent positive Impact 
with the inclusion of satellite data. Kelly et al. (1978) did find a consistent 
positive Impact In a Southern Hemisphere study. This result may have been due 
to the more extensive data poor regions In the Southern Hemisphere. A good sum- 
mary qJ recent model impact studies Is presented by Ohring (1979). 

The Improvement or degradation of numerical model forecasts when satel- 
lite data are combined with conventional data Is often used to imply the quality 
of the satellite data. However, as pointed out by Tracton et al. (1981), th* 
Impact of satellite soundings is a function not only of the satellite data 
quality, but also the analysis and forecast system. Analysis errors may result 
from incompatibility between satellite and conventional data. Consequently, 
analyses based separately on satellite and conventional data may be superior to 
an analysis based on a combined data set. Also, the combination of data seta 
makes the evaluation of the satellite soundings difficult. Analysis techniques 
which take Into account the unique characteristics of the satellite data nay 
produce superior analyses. The model nay also Influence the effect of satellite 
temperature data on numerical model forecasts. If the resolution of the model 
Is too coarse, the model may be Insensitive to small differences in the Initial 
field. Also, inadequate model physics may overwhelm the effects of the changes 
In the Initial field. 

Unlike most previous Impact studies, this study Involves comparisons of 
numerical model forecasts based either entirely on conventional data or entirely 
on satellite data. This approach circumvents problems arising from the mixing 
of data sets. It appears that the only previous st i !y using this approach was 
by Bonner et al. (1976), which produced hemispheric forecasts from analyses 
based only on Nimbus-5 and N0AA-2 temperature soundings. Furthermore, separate 
analyses were produced from TIROS-N and N0AA-6 data allowing intercomparisons 
between two satellite data seta. These analyses were Inserted Into a limited 
area primitive equation model with high horizontal resolution. The resulting 
forecasts are verified against those produced from conventional initial fields. 
From these comparisons, a better understanding of the characteristics and limi- 
tations of satellite data In a forecast situation can be developed. 
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2 • The experiment 

In this study, TIROS-N and NOAA-6 operational satellite soundings were 
used In a numerical forecast model to evaluate the data’s quality and usefulness 
In numerical weather prediction. As noted In the Introduction, this and a com- 
panion study apply the same data seta. In the first, Schmidt et al. (1981) exa- 
mines the satellite sounding quality and the Improvement achieved through manual 
screening of the data. In this second portion of the experiment only the analy- 
ses derived from the 0000 GMT 6 January 1980 TIROS-N and N0AA-6 screened data 
sets are used In a numerical forecast experiment. The data and procedures used 
to construct these satellite-derived analyses are described fully In Schmidt et 
al. (1981). 

To evaluate the model forecast results achieved using TIROS-N and N0AA-6 
analyses, a control forecast was necessary. For this purpose the thickness 
fields obtained from the Limited-area Fine Mesh (LFM) model of the National 
Meteorological Center (NMC) were used. Upper air temperature and geostrophlc 
wind fields were derived using the LFM 1000 mb height analyses and TIROS-N, 

NOAA-6 and LFM thickness analysts. The mean sea level pressure data required by 
the Initialization procedure, were also taken from the 0000 GMT 6 January 1980 
LFM analysis fields. Since the same surface data was used In each case, dif- 
ferences between the satellite and conventional Initial fields were limited to 
the upper air temperature and wind fields. Forecasts out to twenty-four- hours 
were produced from the TIROS-N, NOAA-6 and control initial fields. The control 
run was verified against conventional analyses to evaluate model performance, 
while the forecasts derived from the TIROS-N and NOAA-6 data were compared to 
the control forecasts In order to evaluate the satellite data. This evaluation 
nay provide insights into the ability of the satellite temperature data to 
describe synoptic fields in a numerical forecast situation. It should be empha- 
sized that the satellite data forecasts were not necessarily expected to be 
superior to the control forecasts, particularly since the experiment was con- 
ducted over an area rich In conventional data. Rather, the comparisons were 
Intended to provide Indications of the information content of satellite data. 


3 • Fo recast raodol description 

The Australian Region Primitive Equations model was adapted by the 
Australian Numerical Meteorology Research Centre for use by the University of 
Wisconsin and the National Environmental Satellite Service (NESS). (Hereafter 
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this model will be referred to an the ANMRC model.) The operational version 
upon which the ANMRC model is based is described in McGregor et al. (1978). The 
ANMRC model is a semi-iaplicit primitive equation model with second order trun- 
cation errors in space and time. In the version of the ANMRC model used in this 
experiment, only adiabatic processes were allowed. Additional information on 
this model can be found in Mills et al. (1981). In the vertical, the ANMRC 

model had ten sigma levels while it uses a staggered spatial grid in the hori- 
zontal. 

The ANMRC model grid shown in Figure 1 is on a Lambert coaformal projec- 
tion with a grid spacing of 67.56km at 50° and 20*N. This high resolution not 
only reduces the truncation error, but also is small enough to resolve the 
information contained in the satellite data. Computer memory limitations 
encountered with such a high horizontal resolution restricts the model domain to 
a limited area. Boundary conditions which must then be specified, presented a 
major difficulty in this study. To simplify the experiments, fixed boundary 
conditions vets used. As can be expected, the fixed boundary conditions produced 
large forecast errors near the boundaries. Not only are time variations of the 
atmospheric structure along the boundary neglected by the use of fixed boundary 
conditions, but also an incorrect specification of the initial condition' may 
continuously produce larger errors as the forecast progresses. Thus, only the 
inner region shown in Figure 1 was used to verify 'the forecasts. 

4 « The ANMRC model forcctfsts 

The ANMRC model Initialization procedures required temperature and wind 
analyses at 1000, 850, 700, 500, 300, 250, 200, and lOOnb, moisture data below 
300ab and a sea level pressure field on the ANMRC grid. Since the data were only 
in terns of thickness analyses on the LFM subset grid, the following procedure 
was applied Identically to the screened TIR0S-N, screened NOAA-6 and LFM 
thickness analysis. The thlckneR3 fields derived froa LFM analyses were Included 
to serve as a control forecast In satellite data and model performance evalua- 
tions. The procedure used to produce the required fields was designed to ensure 
consistency between data sets and to minimize the Influence of outside data 
sources . 

First, the lOOOnh height, dew point, temperature and oca level pressure 
analyses were extracted from the LFM analysis and OOhr forecast data sets.'. From 
the TIR0S-N, NOAA-6 - and LFM thickness analyses, temperature fields were derived 













hydrostatically and height fields created using the LFM lOOOob analyses. The sea 
level pressure, temperature, height and dew point temperature analyses were Inter* 
polated to the ANMRC model grid using overlapping quadratic polynomials In two 
dimensions. Since a small area in the southeast corner of the ANMRC model grid Is 
located outside the LFM model grid, an extrapolation to these grid points was 
required. However, in this region any errors introduced by the extrapolation are 
unlikely to significantly effect the forecasts in the verification region. 

On the ANMRC model grid, geostrophic winds verc calculated from the height 
fields. The approximations of geostrophic winds, adiabatic processes, and fixed 
boundary conditions were applied consistently between model runs. Thus, com- 
parisons with the control case should minimize the effects of these approximations 
on the final results. 

The LFM (control), NOAA-6 and TIROS-N data sets produced using the pre- 
ceding procedure were used to initialize and run the ANMRC model on the CRAT-1 
computer at the National Center for Atmospheric Research (NCAR). The model was 
run out to twenty-four hours with a timestep of six minutes. In order to ensure 
consistency between the LFM verifying analyses and the ANMRC model forecasts, the 
filtering in the post processing was altered to simulate NMC's combined filter 
described in Gerrity and Newell (1976). Each forecast was produced consistently 
under the simplified conditions of fixed boundary conditions, geostrophic winds, 
and a minimum of parameterization. Thus, comparisons between ANMRC model runs 
(l.e. the TIROS-N, NOAA-6 and the control run) are more meaningful in data eva- 
luations than comparisons against operational LFM analyses. Model performance 
under these simplified conditions can be evaluated through comparisons between the 
control forecast (the ANMRC model forecast initialized using the LFM analyses) and 
the operational NMC LFM analyses and forecast. 


5. Results 

The TIROS-N, NOAA-6, and control model runs each produced forecasts of four 
different variables at ten levels of the atmosphere for five tine periods (00, 06, 
12, 18, and 24 hours). Since a complete presentation of all these results is 
unman- eable, the evaluation will concentrate on the 00, 12, and 24 hour forecasts 
of the 500nb heights and 700-300ab thicknesses. The 500nb height field was chosen 
for evaluation because of the Importance of this level to forecasters. The 700- 
300mb layer was. chosen for several reasons. Figure 2 shows the effect of the 
ANMRC initialization and post-processing procedures on the.700nb control height 


field. Obviously the errors arising from these procedures are related to the aur- 
f.ace topography, which Ilea below 700ab. The 700-300mb layer la also above and 
below the layera In which the aatelllte temperatures display maximum differences 
from conventional data aourcea. (See Schmidt ct al., 1981.) However, data from 
the lcvela above 700mb and below 300mb may still significantly influence the 700- 
300mb layer forecasts. Also, an inconsistency was found in the model Inltlallra- 
tlon procedure after theae forecasts were run. The result was a large ..eight bias 
in the top sigma level (o - .05) in the forecasts. However, the results suggest 
that the levels below 300mb were not greatly affected. The forecast verification 
will concentrate on the trough Initially located over the Pacific Northwest, since 
the eastern trough is severely affected by the boundary conditions. 


The 500mb level 

The control forecasts are first evaluated against the LFH 500mb analyses to 
give an Indication of the ANKRC model performance under the simplified conditions 
of Initial geostrophlc winds, fixed boundaries and a minimum of parameterization. 
The 00, 12. and 24hr control forecasts along with verifying analyses are shown In 
Figure 3. l n this and subsequent 500mb figures both height and geostrophlc tso- 

tach analyses are shown. 1..e geostrophlc winds are used to indicate the strength 
of the gradient in the height field. 

Since the control forecast was Initialized using the LFH analyses for ‘ 
0000 GMT 6 January, the difference between this analysis (Figure 3a) and the OOhr 
forecast (Figure 3b) indicates the effects of the Initialization and post- . 
processing. A slight decrease ln the areas of geostrophlc wind maxima is ap- 
parent, Indicating a decrease ln the 500mb height gradient. This effect, which la 

least partially due to the filtering done ln the post-processing, should be 
present In nil the raodcl forecasts. 

As the trough over the Pacific Northwest Intensifies and propagates east- 
ward, the 12hr control forecast moves the trough somewhat too quickly and overln- 
tenaif les th, height gradient in the base of the trough. (Compare Figures 3c and 
3d.) However, more apparent differences In the 12hr forecasts are the closed 
center and the nearly north-south trough axis In the control forecast compared to 
the open wave and the positive tilt ln the LFH analysis. As shown In Figures 3e 
and 3f, the location of the trough In the 24hr forecast Is only slightly too far 
to the east. However, the control forecast does not Intensify the trough enough. 

This lack of intensification Is reflected In the weaker geostrophlc jet maximum 
over the central Staten. 
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In order to partially aaaeaa the effect of the almplif lcntlona employed In 
obtaining the control forecast, the KMC LFM forecaata are compared with the control 
lorecasts in Figure 4. Also, since the control Initial fields were taken from the 
LFM analyses, deficiencies in the Initial fields should be reflected In both model 

forecasts* 

Away from the boundaries, two major differences are noted. In the 12hr 
forecasts, the center of the 500mb low over North Dakota Is shifted further south 
In the control forecast, thus Increasing the gradient in the baae of the trough 
(Figurea 4c and 4d). Also, In the 24hr forecasts the LFM model correctly forecasts 
the center of the 500ob low. Compare Figures 4f with 3e. The control forecast 
(Figure 4e) shows a center about sixty meters higher than the LFM forecast. In 
general, the KMC operational LFM forecast and the AHHRC control forecasts were In 
better agreement with each other than either was with the verifying LFM analyses. 
Considering the simplifications used In the control run, it was surprising that the 
control forecasts were In reasonably good agreement wlti the operational LrM 
forecast* 

Since the. satellite forecasts were produced under the same model limitations 
as the control forecasts, many of the same features found In the control forecasts 
can be expected to also be present In the satellite forecasts. Thus, differences 
between the control forecasts and the TIROS-N and NOAA-6 forecasts should reveal 
strengths and limitations of satellite data in a numerical model forecast. The 
companion study by Schmidt ct al. (1981), which gives a more complete description 
of the satellite data and analyses used In the model initialization, found that the 
satellite data tended to be too Warm in troughs and too cold In ridges. Thus, the 
gradients in the height and thickness fields were usually underestimated. In com- 
parisons between TIROS-N and NOAA-6 analyses, Schmidt et nl. found that the riROS-N 
analyses tended to be slightly better than those from NOAA-6. . 

Because of the weaker gradients in the analysis fleld3, one may expect to 
find weaker gradients In the satellite forecasts, with TIROC-N possibly performing 
somewhat better than NOAA-6. The TIROS-N OOhr forecast shows weaker gradients with 
an accompanying underestimation of the geostrophlc Jet maxima in the base of both 
the eastern and western troughs as shown In Figure 5. Also at OOhr tne. western 
trough extended back further Into central Canada. This feature, also founa in the 
NOAA-6 data (Figure 6b), In synoptlcally realistic and nay not have been resolved 
by the widely spaced Canadian radiosonde network. In the 12 and 24hr TIROS-N fore- 
casts the main trough is In approximately the same position -as the control fore- 
cast. However, at 12hrs a weak secondary wave la evident In south-central Canada. 


y 
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Figure 5: Same as Fig. 3 except from the control (a,c,e) and the 

TIROS-N (b,d,f) forecasts. 
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This secondary trough deepens at 24hrs which along with the main trough created a 
very broad trough over the north-central United States. Inspection of the inter- 
mediate forecasts (not shown) reveals that this feature piopagates in 
northern boundary. In both the 12 and 24hr forecasts, the maximum gradients as 
described by the geostrophic jet arc slightly underestimated and too far to the 
south. Also, the overestimation of heights in the center of the 500mb trough seem. 

to propagate through to the 24hr forecast* 

Many features of the NOAA-6 and TIROS-N forecasts are similars the initial 

gradients are too weak, the western trough extends back into Canada in the initial 
data, in the 12 and 24hr forecasts the geostrophic jet maximum was too far to the 
south, and the western trough was too shallow. However, a comparison between 
Figures 5 (TIROS-N) and Figure 6 (NOAA-6) demonstrates some major differences. In 
the initial data, the NOAA-6 analysis shows an anamolous trough located over Mexico 
(Figure 6b). The analysis scheme probably created this anamolous trough by inter- 
polating the gradient created by slightly inferior data into a satellite data gap. 
Apparently, this trough did not greatly influence the forecast as it was nearly 
damped out by 24hr. In the Pacific Northwest, there was little tendency in the 
NOAA-6 forecasts for a secondary trough to form and intensify. (See Figures 6b. d 
and f.) Thus, the troughs in the 12 and 24hr forecasts are considerably more 
realistic than those in the TIROS-N forecasts. The NOAA-6 OOhr fields also 
overestimated the heights in the trough more than the OOhr fields in the TIROS-H 
forecasts. This feature was propagated along with the trough out to 24hrs. How- 
ever, the geostrophic jet maxima were not underestimated in the 12 and 24hr fore- 
casts. Overall, the forecasts using only satellite upper air data did better than 
expected in forecasting the 500mb height field with the NOAA-6 forecasts being 
slightly better in positioning the 500mb trough and the TIROS-N forecasts being 
slightly better in describing the magnitudes of the height fields. 

b. The 700-300mb layer 

The control and verification 700-300mb thickness and thermal wind iootach 
analyses arc shown in Figure 7 and the thickness difference fields between tho cams 
two thickness fields are shown in Figure 8. The effects of the Initialisation on 
the 700-300mb field was once again a reduction of the gradient as indicated by the 
700-300mb thermal wind pattern. The 12hr forecast shows the western thickness 
trough moving too fast resulting in large negative thickness differences over the 
Dakotas and positive differences on the western side of the trough. At both 12 and 
24hrs, the center of ; the maximum thickness gradient is properly located. However. 




Figure 7: 700- 300mh • thickness (solid, dam) and thermal wind isotach 

(dashed, 20 m s~* interval) analyses for the LFM analysis 
(a,e,e) and control forecast (b,d,f). isotach values 
greater than 40 m s ^ are shaded. 
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the 12hr forecast overestimates the gradient at the base of the western trough* At 
24hrs, the thickness trough has become too broad with a slight indication of a 
short wave trough developing on its eastern side* The geos trophic Jet maxlmua 
associated with this pattern has a greater east-west extent* 

Since the TIROS-N and NOAA-6 700-300mb OOhr forecasts shown in Figures 9b 
and 10b have no 40m s~* thermal wind maximum in the western trough, it is apparent 
that the satellite underestimated the thickness gradients* This is also shown in 
Figures 11a and b. The 700-300mb satellite minus control thickness difference 
fields show that for both satellites the western trough axis is initially defined 
by a region of positive thickness anomalies surrounded by negative thickness dif- 
ferences, resulting in the underestimation of the thickness gradient in the 
satellite data. For both TIROS-N and NOAA-6, the resulting overestlmatlon of 
heights are propagated with the trough through 24hr* Despite the underestimation 
of the gradient in the OOhr forecast, the thermal wind jet maxima in the 12 and 
24hr forecasts are of nearly the same magnitude In the satellite forecasts as in 
the control forecast* However, in both the TIROS-N and NOAA-6 forecasts there is a 
tendency for the maximum gradient to be located on the west side of the trough as 
opposed to being located in the base of the control trough. Both satellite fore- 
casts maintain the ridge north of Lake Superior better than the control forecast In 
the 24hr forecasts. (Compare Figures 9f and lOf.) Thus, the trough is not as 
broad as the control forecast yet still broader than the verification. 

Even though both satellite forecasts are similar, certain differences exist* 
In all the TIROS-N forecasts, the magnitude of the thickness field more closely 
resembles the control forecasts and verifying analyses than the NOAA-6 forecast* • 
In the TIROS-N forecast, anomalous troughing Is indicated off Baja, Mexico at 12hrs 
and moves inland at 24hrs. (See Figures 11c and e.) This Is probably due to 
slightly incorrect initial data along the fixed boundary which continuously pro- 
duced errors through the forecast period. At 24hrs both the TIROS-N and NOAA-6 
forecasts have developed indications of a short wave thickness trough over the 
Midwest as did the control forecast. However, the short wave trough in the NOAA-6 
forecast is considerably more pronounced. (Figures 9f and lOf.) 

c. Statistical evaluations 

The statistical evaluation of the forecast was performed using layer mean 
temperatures for the three layers, 1000-700mb, 700-J00mb and 300-100mb. Bias, 
standard deviation and SI scores were calculated. (See Table 1.) The verification 
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Figure 10: Same as Fig. 7 except 

(b,d,f) forecasts. * 
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Table Is Layer mean temperature comparisons between the control, TIROS-N and 
NOAA-6 forecasts and the LFM analyses (top three lines of a, b, and c). Aleo 
compared were the TIROS-N and NOAA-6 forecasts with the control forecasts (last 

two lines of a, b, and c). 

a. BIAS 


1000-7 00MB 


700-300MB 


300-100MB 


00 12 24 00 
HOUR HOUR HOUR HOUR 

CONTROL VS ANAL. .43 .14 -.35 -.36 
TIROS-N VS ANAL. -.15 -.81 -1.83 -.79 
NOAA-6 VS ANAL. .71 .01 -.71 -.90 
TIROS-N VS CONTROL -.58 -.96 -1.48 -.43 
NOAA-6 VS CONTROL .28 -.13 -.36 -.54 


12 

24 

00 

12 

24 

HOUR 

HOUR 

HOUR 

HOUR 

HOUR 

-.03 

-.19 

.26 

-1.69 -2.73 

-1.00 

-1.77 

1.33 

- .81 -1.87 

-.58 

-1.13 

1.44 

- .52 -1.40 

-.97 

-1.57 

1.07 

.88 

.85 

-.54 

- .94 

1.18 

1.18 

1.32 


b. STANDARD DEVIATION 



1000- 

•700MB 


700 

-300MB 


300-100MB 


00 

12 

24 

00 

12 

24 

00 

12 

24 


HOUR 

HOUR 

HOUR 

HOUR 

HOUR 

HOUR 

HOUR 

HOUR 

HOUR 

CONTROL VS ANAL. 

.79 

1.87 

2.22 

.43 

1.33 

1.62 

.50 

1.30 

1.73 

TIROS-N VS ANAL. 

1.44 

2.17 

3.12 

1.17 

1.81 

2.10 

1 .26 

1 *53 

3 .02 

A 7 A 

NOAA-6 VS ANAL. 

2.07 

2.57 

2.81 

1.71 

2.07 

1.74 

1.92 

1.90 

2 #70 

TIROS-N VS CONTROL 

1.32 

1.54 

2.40 

1.03 

1.66 

2.41 

.98 

1 *56 

2.63 

NOAA-6 VS CONTROL 

2.18 

2.37 

2.52 

1.59 

1.87 

1.69 

1.62 

2.02 

2.49 





c. Si 







1000 

-700MB 


700-300MB 



300-IQ0MB 


00 

12 

24 

00 

12 

24 

00 

12 

24 


HOUR 

HOUR 

HOUR 

HOUR 

HOUR 

HOUR 

HOUR 

HOUR 

HOUR 

CONTROL VS ANAL. 

i 

17.0 

33.0 

39.9 

10.3 

28.9 

32.1 

11.7 

31.6 

38.3 

TIROS-N VS ANAL. 

36.0 

43.9 

56.0 

25.0 

34.0 

43.7 

24.7 

34 .8 

56 .2 

NOAA-6 VS ANAL. 

46.2 

51.0 

53.3 

32.4 

36.8 

41.7 

36.7 

40.4 

55.5 

TIROS-N VS CONTROL 

30.2 

30.6 

36.6 

21.6 

34.0 

45.0 

21.0 

37.2 

47 .4 

NOAA-6 VS CONTROL 

41.1 

40.3 

39.4 

28.8 

37.2 

38.7 

31.7 

43.8 

52 .6 


region previously shown in Figure 1 was chosen to be over an area rich in conven- 
tional data and was positioned to minimize the more severe boundary effects. In 
conducting the statistical comparisons , every third grid point was used In order 
to make the results reasonably compatible with those done on the LFK subset grid 
used by Schmidt et al. (1981). Recall that the LFM grid spacing is 190.5km con- 
pared with 67.56km used in the ANMRC model. The satellite forecasts are verified 
against both the LFM analyses and the control run. Verification against the 
control run was intended to reduce the effects on the statistics of the basic 
model errors and consistent errors resulting from the simplified experimental 
design. 

The OOhr comparisons of the control vereus the LFM analyses reveals the 
effect of the initialization and post-processing on the conventional data. For 
all three statistics, the smallest initialization and post-processing errors are 
located in the 700-3G0mb layer. Thus, the eailler choice of the 700-300mb layer 
based on the poor sounding quality below 700mb and above 300mb is further 
justified by the small effects of Initialization and post-processing in this 
layer. 

In the 300-100mb layer during the 24hr forecast period, the bias changes 
by nearly 3°C (from positive to negative values) in all three forecasts when cost- 
pared to the LFM analyses. (See top three lines of Table la.) Closer inspection 
of the data (not shown) reveals that a large part of this bias change was in the 
200-100mb layer. It is likely that much of thi3 feature is due to the previously 
mentioned inconsistency in the initialization of the top sigma level. 

In all cases when the forecasts were compared to the LFM analysis, the 
control forecast had lower standard deviations and SI scores than the satellite 
forecasts. This implies, as expected, that the satellite analyses are less 
accurate than the conventional LFM analysis. However, it is encouraging that the 
standard deviations of the satellite forecasts versus the LFM analyses are oaly 
about 1°C greater than the corresponding control versus LFM analysis values. 

Also, the SI scores for the 12hr forecasts in the 700-300mb layers (lines 2 and 3 
in Table 1c) are certainly respectable. The lowest layer, 1000-700mb, appears to 
have abnormally large standard deviations and SI scores. However, when the 
satellite forecasts are compared to the control forecasts these differences are 
generally reduced. This indicates that much of the difference may be a result of 
the initialization, post-processing or modelling errors, rather than satellite 
data errors. 


-103- 


In comparing the two satellites with each other an interesting feature is 
noted, At 00 and 12hr the TIROS-N forecasts are distinctly superior, but at 24hr, 
the N0AA-6 forecast is somewhat better. Inspection of the 700-300mb thickness 
difference fields In Figure 11 reveals that much of the decrease in the TIROS-N 
forecast skill betweeen 12 and 24 hr Is probably due to the large area of negative 
differences in the Southwest States. As mentioned previously, this error was pro- 
bably due to the incorrect specification of the thermal structure along the fixed 
boundary. In the 12hr forecast the area is still far enough west to be mostly out 
of the verification region. 

In summary, the results involving 500mb height fields, 700-300mb thickness 
fields and statistical parameters are consistent. While the major troughs and 
ridges are reasonably positioned in the satellite forecasts, the troughs are too 
warm and the ridges are too cold. As a result the overall gradients are somewhat 
reduced from those in the LFM analyses and the control forecasts. These patterns 
are similar to those found in the analyses based only on satellite data. (See 
Schmidt et al., 1981.) However, the maximum gradients tend to be of the same 
magnitude as the control forecast, but somewhat misplaced. Also, the forecasts 
unfortunately appear to be influenced by the use of fixed boundaries, especially 
the TIROS-N forecasts. 


6. Conclusions 

Using a high resolution limited area model, twenty-four hour forecasts have 
been produced from initial fields based 'only upon TIROS-N and NOAA— 6 satellite 
upper air data and conventional surface data. These forecasts were evaluated 
against forecasts made using conventional data under the same constraints as used 
In making forecasts based on the satellite data. The evaluation of the Infor- 
mation content of the satellite data in a numerical forecast situation was allowed 
by the complete separation of satellite and conventional upper air data and con- 
sistent preparation of the initial fields. Since two separate initial fields were 
derived using temperature data from two different satellites, an intercomparison 
between the two satellite forecasts was also possible. * ; 

The initial fields as described by the satellite data contained temperature 
gradients which were weaker than those in the conventional analyses. Ab shown by 
Schmidt ct al. (1981) In a more complete study of the basic satellite analyses, 
this situation results from the overeat imat ion of the temperature in troughs .. This 
has been a persistent problem found in many previous studies (e.g. Koehler, 1981, 








-105- 


using Nimbus-6 data and Schlatter, 1981, using TIROS-N data). In the satellite 
forecast, this characteristic of the temperature field was propagated alou^ with 
the synoptic features* However, in the forecasts the maximum temperature gra- 
dients were not generally underestimated but rather they were somewhat incorrectly 
located* Despite the initial underestimation of the temperature gradients,, the 
satellite forecasts were still able to define the location of the major trough and 
ridge posltlono with a reasonable degree of accuracy under fairly severe limita- 
tions* 

The satellite forecasts along with the control forecast were limited by the 
use of initial geoatrophic winds, fixed boundary conditions and very limited pa- 
rameterization. The lack of parameterization is unlikely to have severely 
affected the forecasts. The use of initial geoatrophic winds undoubtedly had a 
detrimental effect on the forecasts* However, since the control forecast was 
similarly degraded, comparisons to the control fo:ecas r . may have reduced the 
detrimental effect of using initial geoatrophic wind.* Unfortunately this is not 
true for the boundaiy conditions* Since fixed boundary conditions were used, an 
initial error in the boundary data may continue to generate errors throughout the 
forecast period. Ah excellent example of this is shown along the western boundary. 
The errors produced by the improper initial boundary values were .'ble to propagate 
well into the Southwest United States in twenty-four hours. Not only were the 
TIROS-N forecasts degraded by the initial data along the weste.*n boundary, but 
also problems probably existed along the northern boundary as shown by the 

appearance of a secondary trough at 500mb. Since LFM, TIROS-N and NOAA-6 analyses 

% 

were available for 1200 GHT 6 January end 0000 GMT 7 January, it nay have been 
better to neat the forecast# using these analyses to update the boundary values. 

Comparisons between the TIROS-N and N0AA-6 forecasts were difficult, since 
it appears that the TIROS-N forecast was more severely degraded by poor initial 
values along the fixed boundary. Statistically the TIROS-N data initially and at 
12hr produced n better forecast than the N0AA-6 data. By 24hr the errors 
resulting from initial errors along the boundar> had propagated Inward far enough 
to severely effect the statistics. Thus, at 24hr the TIROS-N statistics were ef- 
ficiently degraded causing the N0AA-6 forecast statistics to be relatively 
superior. Subjectively, both satellites provided reasonable forecasts considering 
the synoptic situation. The N0AA-6 forecasts generally were slightly better at 
locating the synoptic features, while TIROS-N forecasts better described the 
magnitude of the temperature field. Considering the similarity of the initial 
Zielda it is not surprising that many nspecto of the satellite forecastr were more 
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similar to each other than to the control or verifying analyses* It Is also 
possible that differences in the satellite data and the conventional data vara not 
only a result of satellite data errors , but also the satellite data ray have 
resolved features not found in the conventional data* 

The similarity in forecasts and initial fields further demonstrates that 
the satellite errors are not random but rather highly correlated to the synoptic 
situation* Correlated errors are more likely to be incorporated into the model 
forecasts than random errors* Therefore, the operational satellite data appear 
unlikely to produce a positive Impact on numerical forecasts in a relatively data 
dense region unless the correlated nature of the data errors is either reduced or 
corrected using a thorough understanding of the correlated error otrocture* It 
should be noted that satellite soundings were not intended to completely replace 
conventional data but rather to provide supplementary information in data poor 
regions* In these regions, the satellite data are likely to contain enough addi- 
tional information to improve numerical forecasts, if the difficult problem of 
combining data sources with differing error characteristics can be calved* 
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